PRIMER FOR THE ALGEBRAIC GEOMETRY OF SANDPILES 
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Abstract. The Abelian Sandpile Model (ASM) is a game played on a graph 
realizing the dynamics implicit in the discrete Laplacian matrix of the graph. 
The purpose of this primer is to apply the theory of lattice ideals from al- 
gebraic geometry to the Laplacian matrix, drawing out connections with the 
ASM. An extended summary of the ASM and of the required algebraic ge- 
ometry is provided. New results include a characterization of graphs whose 
Laplacian lattice ideals are complete intersection ideals; a new construction of 
arithmetically Gorenstein ideals; a generalization to directed multigraphs of a 
duality theorem between elements of the sandpile group of a graph and the 
graph's superstable configurations (parking functions); and a characterization 
of the top Betti number of the minimal free resolution of the Laplacian lat- 
tice ideal as the number of elements of the sandpile group of least degree. A 
characterization of all the Betti numbers is conjectured. 



1. Introduction 

This is a primer on the algebraic geometry of sandpiles based on lectures given 
by the first author in an undergraduate Topics in Algebra course at Reed College 
in the fall of 2008 and on subsequent summer and undergraduate thesis projects by 
the second and third authors. It is assumed that the reader has no background in 
algebraic geometry or the theory of sandpiles but is willing to consult introductory 
outside sources such as lOJ and [T7] . 

The Abelian Sandpile Model (ASM) is a game in which one is allowed to stack 
grains of sand on the vertices of a graph G. If a vertex acquires too much sand, 
a grain of sand will be fired to each neighboring vertex. These vertices, in turn, 
may become unstable, and an avalanche of vertex firings may ensue. One vertex 
is usually specified as a sink. Its purpose is to absorb sand fired into it, allowing 
avalanches caused by the addition of sand to eventually come to a halt. The ASM 
associates a group, the sandpile group, to this sand-firing process. The firing rule 
and the sandpile group are intimately connected to the Laplacian of G. 

In algebraic geometry, there is a way of associating a collection of polynomial 
equations to an integer matrix. These polynomials span the lattice ideal corre- 
sponding to the matrix. Our purpose is to apply the theory of lattice ideals in the 
special case where the matrix in question is the Laplacian matrix of a graph G, 
expressing the results in terms of sand on a graph. 

There is another, more widely-known, connection between algebraic geometry 
and sandpiles. It comes from viewing a graph as a discrete version of a Riemann 
surface (i.e., of an algebraic curve over C). As part of this connection, there is a rich 
theory of divisors on graphs, including a version of the Riemann- Roch theorem [3]. 
In Sections 7 and 8, we see that this theory is also relevant for our purposes. 
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We now give a summary of the paper by section. Section 2 is an extended 
outline of algebraic results associated with the Abelian Sandpile Model on a graph. 
What might be new here is a novel treatment of burning configurations (Speer's 
script algorithm), an extension of the result expressing the independence of the 
sandpile group from the choice of sink vertex, and the exposition of the fact that 
an undirected planar graph and its dual have isomorphic sandpile groups. 

After a brief summary of the theory of lattice ideals in Section 3, our main 
object of study — the toppling ideal of a graph — is introduced in Section 4. The first 
paper on the algebraic geometry of sandpiles of which we are aware is Polynomial 
ideals for sandpiles and their Grobner bases, by Cori, Rossin, and Salvy [9 . That 
paper defines the toppling ideal of an undirected graph and computes a Grobner 
basis for the ideal with respect to a certain natural monomial ordering. Sections 4 
and 5 — building on results in the undergraduate thesis of the second author 1 2 7j — 
extend their work, putting it in the context of lattice ideals and, in Theorem |5.11[ 



generalizing the Grobner basis result to the case of directed multigraphs. The proof 
of Proposition 4.2 giving generators for the toppling ideal, is representative of the 
interplay between algebraic geometry and sandpile theory. 

By Theorem |4.11[ any lattice ideal whose zero set is finite is the lattice ideal 
corresponding to some directed multigraph. In that sense, the potential application 
of sandpile methods to lattice ideals is quite broad. As an application of algebraic 
geometry to the ASM, Corollary |5.15 uses Grobner bases to establish a duality 
between elements of the sandpile group and superstable configurations (G-parking 
functions). The result is well-known for undirected graphs. The proof given here is 
the only one of which we know that works in the more general setting of a directed 
multigraph. 

Section 6 gives an explicit description of the zero set of the toppling ideal. It is 
a generic orbit of a faithful representation of the sandpile group of the graph. The 
affine Hilbert function of the toppling ideal is defined in terms of the sandpile group. 
It is related to the Tutte polynomial of the graph by a theorem of Merino pT] . 
Proposition |6.18| shows that the set of zeros of the toppling ideal satisfies the Cayley- 
Bacharach property. 

Section 7 summarizes the Riemann-Roch theory for graphs and includes results 
obtained in the undergraduate thesis of the third author concerning the minimal free 
resolution of the homogeneous toppling ideal of an undirected graph. The resolution 
is graded by the class group of the graph, closely related to the sandpile group. 
By a theorem of Hochster, the Betti numbers are determined by the simplicial 
homology of complexes forming the supports of complete linear systems on the 
graph. By Theorem |7.7[ the top Betti number counts the following structures 
on a graph: the elements of the sandpile group of minimal degree, the maximal 
degree superstable configurations, the maximal G-parking functions, the acyclic 
orientations with a unique fixed source, and the non-special divisors. Conjecture |7.9| 
suggests a characterization all of the Betti numbers in terms of sandpile groups of 
graphs associated with connected partitions (bonds) of the original graph. For 
more on resolutions of toppling ideals and a generalization of the Riemann-Roch 
theory for graphs to certain monomial ideals, see the paper by Manjunath and 
Sturmfels [10]. 

Finally, in Section 8, we characterize directed multigraphs whose homogeneous 
toppling ideals are complete intersection ideals. Further, we give a new method 
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of constructing directed multigraphs whose homogeneous topphng ideals are arith- 
metically Gorenstein. 
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2. Sandpiles 

In this section we summarize the basic theory of sandpile groups. Many results 
are stated without proof. The reader is referred to [T7] and for a thorough 
introduction to the subject. 

2.1. Graph theory. Let G = (V, E) be a directed multigraph with a finite set of 
vertices V and of directed edges E. For e = (u, u) G E' C y x we write :— u 
and := v for the tail and head of e, respectively. If = e^, the edge is a loop. 
These are allowed but do not add much to the theory. By "multigraph" we will 
mean that there is a weight function, 

wt:VxV^N, 

such that wt(M, w) > if and only if (u, v) G E. One may think of an edge e = {u, v) 
of wt(e) as wt(e) edges connecting u to v. For v gV, 



outdeg(w) :— ^ wt(e) 
indeg(v) := ^ wt(e). 



e^E:e+—v 

The graph G is undirected if wt(M, v) — wt(u, u) for all u,v € V, and it is unweighted 
if the weights of all of its edges are 1. If G is undirected, we use the notation 
deg(w) := outdeg(u) = indeg(i;). 

A vertex u is accessible from a vertex v if there is a directed path beginning at u 
and ending at v. A vertex s is globally accessible if it is accessible from all vertices 
of G. Throughout this primer, we v^rill only consider graphs having at least 
one globally accessible vertex. In particular, undirected graphs are assumed to 
be connected. 

Definition 2.1. A sandpile graph is a triple (V, E, s) consisting of a finite, directed 
multigraph {V, E) with a globally accessible vertex s. The vertex s is called the 
sink of the sandpile graph. If, in addition, s has outdegree 0, it is called an absolute 
sink. The nonsink vertices are denoted V := V \ {s}. 
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If G = (y, E, s) is a sandpile graph, we will also refer to the graph (V, E) as G. 
Note that the sink of a sandpile graph need not be absolute; however, for much of 
what we say, one could safely remove outgoing edges from the sink without changing 
the theory. 

Example 2.2. Figure [T] depicts a sandpile graph G. Edges (wi, W2), (w2, wi), (^2, s), 
(w3, s), and (s, ^3) are directed edges with weights 1, 2, 1, 1, 5, respectively; {i^i, W3} 
is an undirected edge of weight 3; and {v2,V3} is an undirected, unweighted edge. 
Although s is the sink of the sandpile graph, outdeg(s) = 5. 




^2 <IF 9f ^3 



Figure 1. Sandpile graph G with sink s. 

For any finite set X, let 

IjX = {X^xex flx 2; : flj. G Z for all x e X} 
be the free Abelian group on X. Restricting to nonnegative coefficients gives NX. 

Notation 2.3. For a,b £ ZX, we define deg(a) — J^xex '^x ^^'^ a > b ii Gx > bx 
for all X (z X . We say a is nonnegative if a > 0. The support of a is 

supp(a) = {x e A" : flj. 7^ 0}. 

Similar notation is used for integer vectors. 

Let G = {V, E, s) be a sandpile graph. 

Definition 2.4. The (full) Laplacian of G is the mapping of groups A : ZV — >■ ZV 
given on vertices v by 

A(t;) := outdeg(u) v — wt(i'. u) u. 

The reduced Laplacian of G is the mapping of groups A : ZV — >■ ZV given on 
nonsink vertices v by 

A(w) := outdeg(w) v — wt(w, u) u, 
uev 

summing this time only over V. 

The Laplacian just defined is dual to the Laplacian one often sees in the litera- 
ture. Define L : ^ by 



L(f){v) := ^ wt{v,u){(j){v) - (j>{u)) 
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for a function (f) G and vertex v. Say V = {vi, . . . ,Vn+i}, and define the 
diagonal matrix D = diag(outdeg(wi), . . . , outdeg(w„+i)). Let A be the adjacency 
matrix, A, given by Aij — wt{vi, Vj). Fixing an ordering ui, . . . , u„+i of the vertices 
identifies with Z""""^ and identifies L with the (n + 1) x (n + 1) matrix 



The matrix for our Laplacian A of G is the transpose of L. 

A spanning tree directed into s is a subgraph T of G with the properties: (1) T 
contains all of the vertices of G, (2) the weight of each edge in T is the same as its 
weight as an edge of G, (3) for each vertex, there is a directed path in T to s, (4) 
for each vertex v ^ s, there is exactly one edge of T whose tail is v, and (5) the 
outdegree of s is 0. If T is a spanning tree directed into s, then its weight, denoted 
wt(T), is the product of the weights of its edges. The following is a basic theorem 
in graph theory. 

Theorem 2.5 (Matrix- Tree) . The determinant of the reduced Laplacian of G is 

the sum, of the weights of all its directed spanning trees into the sink. 

It will occasionally be useful to consider a more restricted class of graphs. 

Definition 2.6. A directed multigraph G = (V, E) is Eulerian if each of its vertices 

is globally accessible and indeg(w) — outdeg(w) for all v <eV . 

Every undirected graph is Eulerian. The condition that indeg(f) = outdeg(i') 
for all vertices v is equivalent to having 1 e ker A. 

2.2. The Sandpile Group. Let G = (V, E, s) be a sandpile graph with nonsink 
vertices V. 

Definition 2.7. A (sandpile) configuration on G is an element of Zy. A configu- 
ration c = ^,,gy c.„v is stable at a vertex w G if < outdeg('u). Otherwise, it is 
unstable. A configuration is stable if it is stable at each v & V . 

As the name suggests, we think of a configuration c as a pile of sand on the 
nonsink vertices of G having Cy grains of sand at vertex v. Sand can be redistributed 
on the graph by vertex firings (or topplings). Firing v & V m configuration c gives 
the new configuration, 



When V fires, we imagine wt(e) grains of sand traveling along each edge e emanating 
from V and being deposited at e+. If e+ = s, then sand sent along e disappears 
down the sink. If c is unstable at v, we say that firing v is legal. The sequence 
of nonsink vertices ui , . . . , Ufe is a legal firing sequence for a configuration c if it is 
legal to fire ui and then it is legal to fire each subsequent Ui from the configuration 
obtained by firing ui, . . . The configuration resulting from applying a legal 

firing sequence to c is the configuration c = c — A a where a € ZV^ is such that (T„ 
is the number of times vertex v appears in the sequence. We write 



L = D-A. 




uev 



= c — Av. 
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In general, we write c — > c if c is the result of applying a legal firing sequence to c. 
In this case, since the reduced Laplacian is invertible (by the Matrix- Tree theorem, 
for instance), there exists a unique a e TiV such that c = c — Act. This a is called 
the firing script or firing vector for c — > c. 

We have the following existence and uniqueness theorem. 

Theorem 2.8. Let c he a sandpile configuration. 

(1) There exists a stable configuration c such that c c. 

(2) Suppose c ^ c with script a and c ^ c! with script a' . Then if c is stable, 
a' > a. If c and c! are both stable, then c = c' . 

Definition 2.9. Let c be a configuration on G. The stabilization of a configura- 
tion c, denoted c° , is the unique stable configuration c such that c — > c. 

Let denote the set of nonnegative stable configurations on G. Then is a 
commutative monoid under stable addition 

a®b:= {a + b)°. 

Thus, stable addition is vector addition in NV followed by stabilization. The iden- 
tity is the zero configuration. 

Definition 2.10. A configuration c is accessible if for each configuration a, there 
exists a configuration b such that a + — > c. A configuration c is recurrent if it is 
nonnegative, accessible, and stable. 

Definition 2.11. The maximal stable configuration on G is the configuration 

Cmax = ^(outdeg(w) - l)t;. 

Proposition 2.12. A configuration c is recurrent if and only if there exists a 
configuration a > such that 

C = a ® Cmax- 

It is not hard to see that the recurrent elements form a semigroup. In fact, they 
form a group. 

Theorem 2.13. The collection of recurrent configurations of G forms a group 
under stable addition. 

Definition 2.14. The group of recurrent configurations of a sandpile graph G is 
called the sandpile group of G and denoted by S{G). 



By Proposition 2.12[ the sandpile group can be found by a systematically adding 
sand to Cmax and stabilizing. Considering a graph consisting of otherwise uncon- 
nected vertices connected into a common sink by edges of various weights, one sees 
that every finite Abelian group is the sandpile group for some graph. 

Example 2.15. The elements of the sandpile group for the sandpile graph in 
Figure [T] are listed below using the notation (ci, C2, C3) :— civi + C2V2 + C3W3: 

(3,3,4) (3,3,3) (3,2,4) (2,3,4) (3,3,2) (3,2,3) (2,3,3) 
(3,1,4) (2,2,4) (1,3,4) (3,2,2) (2,2,3) (1,3,3) (3,0,4) 
(2,1,4) (1,2,4) (0,3,4) (1,2,3) (0,3,3) (2,0,4) (1,1,4) 
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Although the zero configuration is the identity for 7W, it is seldom the identity 
for S{G). The following is an easy exercise. 

Proposition 2.16. The following are equivalent: 

(1) the zero- configuration is recurrent; 

(2) every stable configuration is recurrent; 

(3) every directed cycle of G passes through the sink vertex. 

We now give another description of the sandpile group. 

Definition 2.17. The Laplacian lattice, C C ZV, is the image of A. The reduced 
Laplacian lattice, £ C ZV, is the image of A. The critical group for G is 

C(G) = ZV/C. 

Theorem 2.18. There is an isomorphism of Abelian groups 

S{G) C{G) 
c I— >■ c + C. 

Thus, each element of ZV is equivalent to a unique recurrent element modulo 
the reduced Laplacian lattice. The identity of the sandpile group is the recurrent 
configuration in £. It can be calculated as 

^ — ((Cmax (2Cniax) ) ^~ Cmax) ■ 



Note that rj = mod £, and since Cmax— (2c„iax)° > 0, Proposition 2.12 guarantees 
that 77 is recurrent. 

Example 2.19. The reduced Laplacian of the sandpile graph in Figure [l] is 

/ 4 -2 -3 
A = -1 4-1 

\ -3 -1 5 

The Smith normal form of A is diag(l, 1, 21). Hence, S{G) « Z/21Z. The identity 
is (3,1,4), computed as follows: 

(C,„ax - (2c„,ax)°) + Cn,ax = ((3, 3, 4) - (6, 6, 8)°) + (3, 3, 4) 

= ((3, 3, 4) -(2, 0,4)) + (3, 3, 4) 
= (4,6,4)^(3,1,4). 

As a consequence of the Matrix- Tree theorem, we have the following. 

Corollary 2.20. The order of S{G) is the sum of the weights of G's directed 
spanning trees into s. 

Remark 2.21. Babai [T] has noted another characterization of the sandpile group: 
it is the principal semi-ideal in M generated by Cmax, which turns out to be the 
intersection of all the semi-ideals of A4 . 

Remark 2.22. In the literature, a sandpile configuration is often taken to be an 
element of Z^. We prefer to work in the dual group ZV = Hom(Z^,Z) so that 
the functor that takes a sandpile graph to its sandpile group is covariant. Suppose 
that G = {V, E, s) and G' = (V, E' , s') are sandpile graphs with reduced Laplacian 
lattices C and respectively. Let 5" : G' — > G be a mapping of graphs that 
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maps s' to s. Applying homz( • ,Z) to the natural induced map — > yields 
; ZV' ZV. If ^(/C') C £, there is an induced mapping of sandpile groups. 
This condition would seem to define a reasonable set of morphisms, then, for a 
category of sandpile groups. For work on the category theory of sandpile groups, 
see [5] and |31j . For the notion of a harmonic morphism of graphs, see [1]. 

2.3. Superstables. Let c = u + u be a configuration on the (unweighted, undi- 
rected) sandpile graph in Figure ^ with sink s. The vertices u and v are both 




s 



Figure 2. Graph G. 

stable in c, so there are no legal vertex firings: firing cither vertex would result in a 
negative amount of sand on a vertex. However, firing both vertices simultaneously 
results in a nonnegative configuration, the zero configuration. Each nonsink vertex 
loses two grains of sand, but each also gains a grain from the other. 

Definition 2.23. Let c be a configuration on the sandpile graph G = (V, E,s). A 
script-firing, also called a cluster- or multiset- firing, with (firing) script a S NV is 
the operation that replaces c with c — Aa. The script-firing is legal if cr > and 
(c — Acr),u > for each v £ swpp{a). Thus, if c > 0, the script-firing with script 
(T > is legal if and only if c — A cr > 0. 

A configuration c is superstable if c is nonnegative and has no legal script-firings. 

The idea of a G-parking function is essentially the same as that of a superstable 
configuration: 

Definition 2.24. Let G = (V, E, s) be a sandpile graph. A G-parking function [2H] 
(with respect to s) is a function f : V Z such that there exists a superstable 
configuration c on G with the property that f{v) — Cy for v € V and /(s) ~ —1. 

An acyclic orientation of an undirected graph G is a choice of orientation for 
each edge of G such that the resulting directed graph has no directed cycles. A 
vertex w is a source for an acyclic orientation if all the edges incident on v are 
directed away from v. If O is an acyclic orientation and v € V, then indegQ(u) 
denotes the indegree of v for the directed graph corresponding to O. 

Theorem 2.25 ([5]). Let G = (V, s) be an undirected sandpile graph. Then 
there is a bijection between the set of acyclic orientations of G with unique source 
s and the set of superstable configurations of G of highest degree. If O is an acyclic 
orientation, the corresponding maximal superstable configuration is given by 

^(indegc,(D) - 1) V. 

For an extension of the previous theorem from maximal superstable configura- 
tions to all superstable configurations (and a connection with hyperplane arrage- 
ments), see [T5] . 
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2.4. Burning configurations. Speer's script algorithm [29] generalizes the burn- 
ing algorithm of Dhar, testing whether a configuration is recurrent. We present a 
variation on Speer's algorithm using burning configurations. 

Definition 2.26. A configuration 6 is a burning configuration if it has the following 
three properties: 

(1) b e C, 

(2) 6 > 0, 

(3) for all V E V, there exists a path to v from some element of supp(6). 

If 5 is a burning configuration, we call ab — (A)^^6 the script or the firing vector 
for b. 

Theorem 2.27 f|26j). Let b be the burning configuration with script df,. 

(1) {kb)° is the identity configuration for k ^ 0. 

(2) A configuration c is recurrent if and only if the stabilization of c + b is c. 

(3) A configuration c is recurrent if and only if the firing vector for the stabi- 
lization of b + c is ab- 

(4) at > 1. 

(5) If c is a configuration and t is the firing vector for the stabilization ofc + b, 
then T < ab- 

Thus, a configuration c is in the sandpile group if and only if adding a burning 
configuration to c and stabilizing returns c, or if, equivalently, the firing script for 
the stabilization is equal to the burning script. For the case of an undirected graph, 
as we see in the following theorem, one may take 1 as the firing script. Adding the 
burning configuration to a configuration c in that case can be thought of as placing c 
on the graph, then firing the sink vertex. Checking whether each vertex fires exactly 
once in the subsequent stabilization is known as Dhar's algorithm. 

Theorem 2.28 (|29).pS]). There exists a unique burning configuration b with script 
(Tfc = A^^6 having the following property: if ab' is the script for a burning configu- 
ration b' , then ab' > ab- For this b, we have: 

(1) For all V V , b^ < outdeg(w) unless v is a source, i.e., unless indeg(u) = 0, 
in which case by — outdeg(D). Thus, b is stable unless G has a source, and 
in any case, by < outdeg(u) for all v. 

(2) CTfe > 1 with equality if and only if G has no "selfish" vertices, i.e., no vertex 
v eV with indeg(w) > outdeg(w). 

We call this b the minimal burning configuration and its script, ab, the minimal 
burning script. 

Remark 2.29. To compute the minimal burning configuration, start with b equal to 
the sum of the columns of A. If > 0, stop. Otherwise, if 6^ < for some v G V, 
replace b hy b + A{v). Repeat until 6 > 0. 

Example 2.30. We would like to compute the minimal burning configuration and 
corresponding script for the sandpile graph G in Figurejl] Continuing Example 2.19 



the sum of the columns of A is (—1,2,1)*. Since the first entry of the sum is 
negative, add in the first column of A to get (3, 1, —2)*. Since the third entry is now 
negative, add in the third column of A to get (0, 0, 3). Thus, the minimal burning 
configuration is b = (0,0,3), and the burning script is ab = (2,1,2), recording the 
columns of A used to obtain b. 
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2.5. Some isomorphisms. 

2.5.1. Choice of sink vertex. Lemma 4.12 of [171 states that for Eulerian graphs, 
the sandpile group is, up to isomorphism, independent of the choice of sink. Here, 
we present a generahzation of that result. 

Let G = (y, E, s) be a sandpile graph. Recall that C{G) :— ZV/C is the critical 
group of G, isomorphic to the sandpile group, S{G), by Theorem 2.18 Let 



ZVo := {ceZV : deg(c) = 0}. 

Since the image of the Laplacian A is contained in Zuq, we may define the mapping 
Ao-.ZV ^ ZVo by Ao(c) A(c) for aU c e ZV. 

Proposition 2.31 ( 26 ). 

(1) There is a commutative diagram with exact rows 



> ZV > ZV - 



->0 



> ZV/ ker A ZVq > C{G) 



^0. 



where l{v) := v + ker A and e — v — s for all v £ V . 

(2) For each v £ V , let Ty be the sum of the weights of all spanning trees 
directed into v, let d = gcd{r„ : u S V} , and let fy := Ty/d. Define 
f := J2y € ZV. Then 

ker A = Span^jf}. 

(3) There is a short exact sequence 

— > Z/i,Z — > C{G) — > ZVo/C — > 0. 

Corollary 2.32. If G is an Eulerian graph (in particular, if G is undirected), then 
the sandpile group for G is independent of the choice of sink vertex. 

Proof. Suppose G is Eulerian. Then each vertex is globally accessible. So it makes 
sense to talk about the sandpile group of G with respect to any of its vertices. 
Since indeg(w) = outdeg(w) for all v G V, we have that 1 G ker A. It follows from 
Proposition 2.31| ^ that fy — 1 for all v. Fix a vertex s and consider the sandpile 
group of G with respect to s. It is isomorphic to the critical group (with respect to 
s), and hence isomorphic to ZVq/C by Proposition 2.31| ^ 
not depend on the choice of a sink. 



However, ZV^j L does 
□ 



2.5.2. Planar duality. Let G = {V, E) be an undirected graph. Fix an orientation O 
of the edges of G. Thus, for each {w, u} £ E we have that either {u,v) or {v,u) 
is in O, but not both. Let e — {u,v} £ E, and suppose that {u,v) £ O. In 
the free abelian group ZE, we identify (m, v) with e and (v, u) with — e. We also 
define e~ :— u and e+ := v. 

The (integral) cycle space, C — Cq C ZE, is the Z-span of the cycles of G. 

Example 2.33. Let G be the (undirected) triangle with edges oriented as in Fig- 
ure [sj The cycle space for G is the Z-span of the cycle {x, y) + (y, z) — (a;, z). 
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Figure 3. A triangle with oriented edges. 

For each U C V, define the corresponding cut-set, c^, to be the collection of 
edges of G having one endpoint in U and the other in the complement W^. For 
each e € E, define the sign of e in a cut-set by 

( 1 if e- e t/ and e+ e C/^ 
cr(e, cIj) := I 1 if e- e U" and e+ G U, 
[ otherwise. 

We then write = J2eeE '^(^^ '^u) ^ ^ "^^^ Z-span of the cut-sets of G is the 
(integral) cut space for G, denoted C* . If J7 = {z;} for some v G V, then c* := is 
called a vertex cut. It is well-known that the vertex cuts form a Z-basis for C* . 
Define the boundary mapping by 

e I— > — e~. 

for e G E. We have the following well-known exact sequence (recalling that we are 
assuming G is connected): 

^ C > ZE ZV Z ^ 0. 

A straightforward calculation shows that for each v E V, 

diO = A{v). 

We have the following theorem. 

Theorem 2.34 ([!]). Let G be an undirected sandpile graph. Then 

S{G)kZE/{C + C*). 
The following result appears in 8 . 

Corollary 2.35. Let G be an undirected planar graph, and let G* be its dual. 

Choosing any vertices to serve as sinks, there is an isomorphism of sandpile groups 

S{G)kS{G*). 

Proof. An orientation of G induces a dual orientation on G*: if F and F' are 
adjacent faces in G (vertices of G*) intersecting along edge e, we orient the edge 
e* := F'} of G* as (F, F') if F is to the right of e as one travels from e~ 
to e+ . Sending e to e* then defines an isomorphism ZE ZE* where E* denotes 
the edges of G* . It is well-known that under this isomorphism the cycle space 
(resp., cut space) of G is sent to the cut space (resp., cycle space) of G* . The 
result then follows from Theorem |2.34[ The choice of sink vertices is irrelevant by 
Proposition [23T| □ 
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Remark 2.36. 

(1) The independence of the sandpile group of G, up to isomorphism, of the 
choice of sink is also a consequence of Theorem |2.34[ 

(2) Theorem [2.34 suggests a definition of the sandpile group for an arbitrary 
matroid ([19,). 

(3) As noted in if two undirected (connected) graphs are 2-isomorphic, then 
their corresponding matroids are isomorphic. (See [24j for the definition of 
2-isomorphism and a proof of the Whitney's 2- isomorphism theorem.) It 
then follows from Theorem |2.34| that the sandpile groups for the two graphs 
(having chosen sinks) are isomorphic. 

3. Lattice ideals 

Our reference for this section is [22]. Let A be a finitely generated Abelian 
group, and let ai, . . . ,a„ be a collection of elements generating A. Let Q be the 
subsemigroup of A generated by ai, . . . , a„. In the case where A is finite — the case 
of special interest to us — we have that Q = A. Define ^: Z" — >■ Q by (f>{ei) — Ui, 
and denote its kernel by A. Let {ta ■ a G Q} be indeterminates, and let 

C[Q] = Spanc{ia : a e Q} 

be the group algebra of Q; hence, tatb = ta+b for elements a, 6 e Q. Letting 
R :— C[a;i, . . . , a;„], define a surjection of rings 

^: R ^ C[Q] 

For c € N", we define = Hi^?- Then ip{x'^) is the group algebra element tf,, 
where b = X^f^i c^Oi. 

For u £ Z", we write u = — u~ with , u~ € N" having disjoint support. 

Theorem 3.1. 

(1) The kernel of ijj is the lattice ideal 

/(A) Spanc{x" - : m, i; e N", u - i; G A}. 

(The vector space span, above, forms an ideal.) Hence, tp induces an iso- 
morphism of C-algebras, R/I{A) « C[(3]. 

(2) If £i, . . . , £k are generators for the Z-module, A, then /(A) is the saturation 
of 

J= (x^i" -x^' ■.i^l,...,k) 
with respect to the ideal generated by the product of the indeterminates, 
]Xi=i^i- Thus, 

/(A) ^{f eR: (nr=i ^^)"'f e J for some m e N}. 

(3) The Krull dimension of R/I{A) is n — dim^ A. 

Let U C N" such that X :— {x" : u G U} is a, C-vector space basis for R/I{A). 
Letting g :— (ai, . . . ,a„) G A", 

V^(A) = {tu.g -.ueUj^ita-.aeQ}, 

the last equality holding since R/I[A) and C[(3] are isomorphic as vector spaces 
via ip. Now assume that ^ is a finite group, so that Q = A. Then, ij) induces a 
bijection of X with A, which endows X with the structure of a group isomorphic 
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to A. For u,v £ U, we define x^x" = x"^ where w is the unique element of U for 
which w ■ g — {u v) ■ g. 

A choice of a monomial ordering on R gives a natural choice for [/, namely, those 
u € N" such that is not divisible by the initial term of any element of /(A), e.g., 
not divisible by the initial term of any element of a Grobner basis for /(A). This 
will be discussed in S|5| 

Example 3.2. Let A = 1/21 x Z/3Z with generators Oi = (1,0), = (0,1), 
and aa = (1, 1). The kernel A of ^ A is spanned by (2,0,0), (0,3,0), and 

(1, 1, —1). Hence, the saturation of the ideal {x\ — 1, — 1, xiX2 — x-^ gives the 
lattice ideal /(A). Using a computer algebra system, one computes 

/(A) = {x\ - 1, Xi_X2 - xz, xxxz -x-2,x\- x\,X2x\ -\,x\- Xi). 

By Theorem |3.1 1 ([S]), one expects a finite set of solutions over C to the equations 
formed by setting the generators of /(A) equal to zero — there are six. One vector- 
space basis for R/I{K) is 

l,Xi,X2,X'i,X2Xs„xl. 

4. Toppling ideals 

Let G be a sandpile graph. Identify its vertices with {l,...,n-|-l}, where n + 1 
represents the sink. To avoid ambiguity, we will sometimes denote vertex i by w^. 
By ordering the vertices, we thus have the exact sequence for the sandpile group 
of G, 

^ Z" A Z" ^ 5(G) ^ 0. 
Recall our notation for the reduced Laplacian lattice: 

C = im(A) = kcr(Z" ^ 5(G)). 
Definition 4.1. The toppling ideal for G is the lattice ideal for £, 

/(G) := Spanc{x" — : u — v mod C} d R — C[a;i, . . . ,a;„]. 
The coordinate ring for G is R/ I{G). 

Thus, by Theorem [sl] , we have the isomorphism of C-algebras: 

R/I{G) w C[5(G)]. 
For each nonsink vertex i, define the toppling polynomial 

outdcg(z)— wt(i,z) -p-r wt(i,j) 



Proposition 4.2. The ideal I {G) is generated by the toppling polynomials, {ti\^^i, 
and the polynomial a;'' — 1 where b is any burning configuration. 

Proof. Let J = {t^ : i = 1, . . . ,n) + {x^ — 1). It is clear that J C /(G), and by 
Theorem 3.1| (2|), /(G) is the saturation of J with respect to the ideal (a;i • • • x„). So 



it suffices to show that J is already saturated with respect to that ideal. Suppose 
that {xi ■ ■ ■ Xn)''f S J for some f E R and for some k. For each positive integer to, 
consider the monomial x™^. We think of this monomial as a configuration of sand 
with mbi grains of sand on vertex i. If vertex i of this configuration is unstable, 
we think of firing the vertex as replacing a;™''* by x™'''"'*' rij^^i xj^^^'''\ Performing 
this replacement in x™'' gives an equivalent monomial modulo J. Recall that every 
vertex of G is connected by a directed path from a vertex in the support of b. Thus, 
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by taking m large enough and firing appropriate vertices, we arrive at a monomial 
^ equivalent to x"^^ modulo J and corresponding to a configuration with at least k 
grains of sand at each vertex. Write x'^ = x^{xi ■ ■ for some monomial x^ . 

Modulo J, we have 

= {xi---xnff 

= X^{xi---Xnfl 

= 

= a;"V 

- /• 

Thus, / e J, as required. □ 

Remark 4.3. As in the proof of the above theorem, we can identify a monomial x°' 
with the configuration a on G. If a — >■ 6 as sandpile configurations, then = x^ in 
R/I{G). 

Remark 4.4. The toppling ideal was introduced by Cori, Rossin, and Salvy [S]. 
They considered only undirected graphs and defined the ideal via generators. For 
an undirected graph, the all-Is vector is a burning script, so Proposition |4.2| shows 
that our definition coincides with theirs in the case of an undirected graph. 




Figure 4. Sandpile graph G with sink v^. 



Example 4.5. The sandpile graph G in Figure|4]has a burning script a = (1, 2, 1) 
and corresponding burning configuration b — (0, 1, 2). Thus, 

I[G) — {x\ — X2Xz, x\ — xi,x% — x\,X2x\ — 1). 

Definition 4.6. Let f E R ^ C[xi, . . . ,n], and let Xn+i be another indeterminate. 
The homogenization of / with respect to Xn+i is the homogeneous polynomial 

If / C i? is an ideal, the homogenization of / with respect to Xn+i is the ideal 

l'^ -.^ if'^-.fe I). 

Now consider the exact sequence corresponding to the full Laplacian, 

z"+i A z"+i ^ r'+^/c 

recalling the notation for the Laplacian lattice, £ := im(A). Let S = C[xi, . . . , Xn+i] 
and consider the lattice ideal for C. We here introduce the homogeneous version of 
the toppling ideal. 
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Definition 4.7. The homogeneous toppling ideal for G is 

Ih{G) :— Span£{x" — x" : u — v mod £} C S* = C[xi, . . . , a;„+i]. 

The homogeneous coordinate ring for G is S/Ih{G). 

The following proposition is straightforward. Its hypothesis is satisfied for any 
Eulerian graph and, in particular, for any undirected graph. Moreover, given any 
sandpile graph with sink s, removing all out-edges from s creates a new sandpile 
graph with the same sandpile group and for which the hypothesis of the proposition 
holds. 

Proposition 4.8. // A(i;„+i) e Span2{A(wi), . . . , A(w„)}, then Ih{G) = I[Gf. 

Example 4.9. The graph G in Figure [5] does not satisfy the hypothesis of Propo- 
sition |4.8[ Regarded as a sandpile graph with sink vi , the toppling ideal for G is 
{x\ — 1). As a sandpile graph with sink W2, its toppling ideal is — !)• Its homo- 
geneous toppling ideal is Ih{G) = (xi — X2), equivalent to that of the undirected 
graph with a single edge connecting vi and V2 (or equivalent to that of the directed 
graph consisting of a single directed edge connecting vi to V2)- 



Figure 5. Graph G. 



Remark 4.10. In general, homogenizing the generators of an ideal does not pro- 
duce a complete set of generators for the homogenized ideal. For instance, the 



graph in Example 7.8 has toppling ideal generated by 4 polynomials, whereas its 



homogeneous toppling ideal is minimally generated by 6 polynomials. 

Theorem 4.11. Let C be any submodule ofZ^^ having rank n. Then there exists a 
sandpile graph whose reduced Laplacian lattice is C Every lattice ideal defining a 
finite set of points is the lattice ideal associated with the reduced Laplacian of some 
sandpile graph. 



Proof. In light of Theorem 3.1 (3]), it suffices to prove that given an n x n matrix M 



of rank n, there exists a matrix M' with the same integer column span as M 
and which is the reduced Laplacian matrix of some sandpile graph. Recall that 
a matrix M' is the reduced Laplacian of a directed multigraph if and only if (i) 
deg(c) > for each column c of M', (ii) Af^ > 0, (iii) < for z 7^ j. (If c is a 
column vector of a matrix, then deg(c) is the sum of the entries of c.) If in addition 
M' has full rank, then its corresponding graph has a globally accessible vertex by 



the Matrix- Tree Theorem. The desired matrix M' is produced by Algorithm 4.12 
stated below. It proceeds in three steps, modifying the columns of M using only 
invertible integral column operations. 

First, since M has rank n, not all columns have deg(c) — 0. Using the Euclidean 
algorithm, by adding multiples of one column to another, we set deg(c) to for all 
but one column c of M (line 1). By possibly moving and negating that column, we 
have that deg(ci) = for all but the first column Ci, for which deg(ci) > 0. 
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Next, we repeat the Euclidean algorithm another (n— 2) times, now on the super- 
diagonal entries of each of the first (n — 2) rows in turn (lines 2-9). Again by adding 
multiples of one column to another, we have every entry more than one row above 
the diagonal set to 0. Note that since this step only involves addition of columns 
whose degree is already zero, the column degrees are not affected. Additionally, 
since M had rank n and the last (n — 1) columns have degree zero, we have that 
each of these columns has a nonzero superdiagonal entry. Now by negating columns 
where necessary, we may assume that the nonzero superdiagonal entry of each 
column is negative. 

At this point, the last column satisfies (i) (iii). Assuming the last r columns 
Cn-r+i, ■ ■ ■ ,Cn Satisfy (i)-(iii) for r < n — 2, wc c;laim that for any 1 < ,s < r there is 
a vector e Span2{c„_r+i, . . • , c„} with < and v^_g > and with all other 
entries zero. For r = 1, the vector is obtained by negating c„, so we proceed 
by induction on r. With the hypotheses satisfied for some r, we already have 
appropriate vectors v^,. . . ,v'^~^. To obtain v'^, note that — c„_r+i has a positive 
entry in row (n — r), so by adding appropriate multiples of the for s < r, we 
produce the desired column vector. 

Given that such vectors exist, it is clear that we may iteratively correct the 
columns from right to left by adding multiples of the higher indexed columns. We 
now give this algorithm explicitly. In what follows, v[j] denotes the j-th entry of 
the column vector v, and the Euclidean algorithm terminates when run in-place on 
some set of integers, S, once a single element of S equals the positive GCD of the 
elements of S and every other element of S is zero. 

Algorithm 4.12. 

Input: An n x n matrix M of rank n with columns Ci, . . . , c„. 
Output: The reduced Laplacian matrix A{G) of a directed multigraph G such 
that A{G) = MU for some invcrtiblc integral matrix U. 

1 Run the Euclidean algorithm on the set S = {deg(cfe)} by subtracting one 
column from another at each step. Swap columns so that deg(ci) = gcd(S') and 
deg(ci) = for I > 1. 

2 for k ^ 2 up to n — 1 do 

3 Run the Euclidean algorithm on the set S = {ci[k—\\ : i> k} by subtracting 

one column from another at each step. 

4 Swap columns so that Ck[k — 1] = gcd(S') and Ci[A; — 1] = for i > fc 

5 Cki Cfc 

6 end for 

7 if c„[n — 1] > then 
9 end if 

10 for k ^ n — 1 down to 1 do 

11 for i fc + 2 up to n do // this loop is not entered until k <n — 2 

12 while Cfe [i — 1] > do 

13 Cfc ^ Cfe + Ci 

14 end while 

15 end for 

16 V i Cfc+i 

17 for i ^ k + 2 up to n do // this loop is not entered until k <n — 2 
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18 V ^ \ct[i — l\ \ ■ V + v[i — IJ • d 

19 end for 

20 while Ck[k] < or Ck[n] > do 

21 Cfc ^ Ck+V 

22 end while 

23 end for 

24 return [ci • • • c„] 



□ 



For the sake of the following corollary, a weighted path graph P — ui . . .Uk is a 
graph with vertex set {ui, . . . , Wfe} and weighted edges {(wi, Wi+i) : 1 < i < A:}. If F 
and F' are weighted digraphs, their graph sum is the graph F + F' whose weighted 
adjacency matrix is the sum of those for F and F' . 

Corollary 4.13. Let G be a sandpile graph with vertex set V = {wi, . . . , 
and sink w„+i. Then there exists a weighted path graph P — ti„w„_i • • • viVn+i o,nd 
a directed acyclic graph D on the nonsink vertices V oriented from lower-indexed 
vertices to higher such that the graph sum G' = P + D has the same Laplacian 
lattice as G. 

The above simply states the form of the graph given by the output of Algo- 
rithm 4.12 The graph G" of Corollary 4.13 is not uniquely determined. For in- 



stance, by iterating line 21 of Algorithm 4.12 more times than necessary, one may 
generate infinitely many graphs G" of the form described in the corollary, each with 
Laplacian lattice £. 



Example 4.14. One sandpile graph of the form given by Corollary 4.13 with the 



same Laplacian lattice as the sandpile graph G from Example 4.5 is G' appearing 
in Figure [6| 




V3 




1 ►« 



V2 



Vl 



Figure 6. The sandpile graph G' for Example 4.14 



Question 4.15. When is it the case that a submodule of Z" with rank n is the 
reduced Laplacian lattice of an undirected graph? It is not always the case. For 
instance. Figure |8] is a directed sandpile graph whose lattice ideal is Gorenstein 



(cf. |8]) and with sandpile group of order 5. By Theorem 8.28 any undirected 
graph with Gorenstein lattice ideal must be a tree and would thus have sandpile 
group of order 1. 

5. GrOBNER BASES OF TOPPLING IDEALS 

We recommend [TU] as a general reference for the theory of Grobner bases needed 
in this section. Let i? = C[a;i, . . . , a;„]. 
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Definition 5.1. A monomial order, >, on i? is a total ordering on the monomials 
of R satisfying 

(1) If > x^ then a;=+» > x'=+'' for all c > 0; 

(2) 1 = is the smallest monomial. 

Example 5.2. The following are the most common examples of monomial orders: 

(1) Lexicographic ordering, lex, is defined by > if the left-most nonzero 
entry of a — 6 is positive (i.e., more of the earlier indetcrminatcs). 

(2) Degree lexicographic ordering, deglex, is defined by > x > 6 if deg(a) > 
deg{b) or if deg(a) = deg(6) and the left-most nonzero entry of a — 6 is 
positive (i.e., order by degree and break ties with lex). 

(3) Degree reverse lexicographic ordering, grevlex, is defined by x"* > x'' if 
deg(a) > deg{b) or if deg(a) = deg(6) and the right-most nonzero entry of 
a — 6 is negative (i.e., order by degree then break ties by checking which 
monomial has fewer of the later indetcrminates). 

A monomial multiplied by a constant is called a term. Once a monomial ordering 
is fixed, write ax"' > /3x^ for two terms if a and /3 are nonzero and x" > x^. Each 
/ e i? is a sum of terms corresponding to distinct monomials. We denote the leading 
term — the largest term with respect to the chosen monomial ordering by LT(/). 

Definition 5.3. Fix a monomial ordering on R and let f,gGR. The S -polynomial 
for the pair (/, g) is 

^_ lcm(LT(/),LT(g)) ^ lcm(LT(/), LT(g)) 

Definition 5.4. Fix a monomial ordering on R, and let / be an ideal of R. A finite 
subset r of / is a Grobner basis for / with respect to the given monomial ordering 
if for all / e / there is a ,g e F such that LT(,g) divides LT(/). 

Let F = {gi, ■ ■ ■ ,gm} be the Grobner basis for an ideal ICR with respect to 
some monomial ordering, and let / G M. If / has a term m divisible by LT{gi) for 
some i, then replace / by / — jj^j^ gi- A standard result in the theory of Grobner 
bases is that by repeating this process one arrives at a remainder r that is unique 
with respect to the property that (i) r = f + g for some g € I and (ii) r has no 
terms divisible by any leading term of an element of F. We call this remainder the 
reduction or normal form of / with respect to the Grobner basis F. 

Notation 5.5. The reduction of / with respect to F is denoted by / % F. \i g & R, 
we write f%g for the special case in which / = [g) and F = {g}. 

Proposition 5.6. Fix a monomial ordering on R, and let I be an ideal of R. The 

following are equivalent for a finite subset F of I : 

(1) F is a Grobner basis with respect to the given ordering; 

(2) there is an equality of ideals: (LT(g) : g gT) = (LT(/) : / e I); 

(3) each f £ I may be reduced to by T, i.e., f %T — 0; 

(4) for all g,g' € T, the S -polynomial S{g,g') reduces to by T and T is a 
generating set for I. 

The last criterion is essentially Buchberger's algorithm for calculating a Grobner 

basis: start with any generating set for /, and if / := S{g, g') % F 7^ for some pair 
of generators g and <?', add / to the set of generators and check the S'-pairs again. 
The process eventually stops. 
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Definition 5.7. Fix a monomial ordering on R and let / be an ideal of R. The 
set of monomials of R that are not divisible by the leading term of a Grobner basis 
element for / with respect to the given ordering is called the normal basis for R/I. 

By Macaulay's theorem (Theorem 15.3, [T3]), a normal basis is a vector space basis 
for R/I. 

We now introduce an appropriate monomial ordering for sandpiles, due to Cori, 
Rossin, and Salvy, [S]. 

Definition 5.8. Let G be a sandpile graph with vertices {vi, . . . ,Vn+i} and with 
sink Vn+i- A sandpile monomial ordering on R ~ C[xi, . . . ,Xn] is any grevlex 
ordering for which Xi > xj if the length of the shortest path from vertex Vj to the 
sink is no greater than that for Vi. Given a sandpile monomial ordering > on R, 
the sandpile monomial ordering on S* = C[a;i, . . . ,a;„_|_i] compatible with > is the 
grevlex order extending > for which x^ > Xn+i for i — 1, . . . ,n. 

Proposition 5.9. With notation as in Definition \5.8\ let > be a sandpile monomial 
ordering on R, extended to a compatible sandpile monomial ordering on S. Let 
I <Z R be the toppling ideal for G. 

(1) Let r a Grobner basis for I with respect to >, and let be the subset of S 
formed by homogenizing each element ofT. Then V'' is a Grobner basis for 
the homogenization C S*. 

(2) The normal bases for R/I and for S/ + (xn+i)) consist of the same set 
of monomials. Hence, R/I and S/{I^ + (a;„+i)) are isomorphic as vector 
spaces. 

Proof. The first part of the proposition is a general result for grevlex orderings 
(cf. Exercise 5, §8.4, jTU]). It is straightforward to check that if / G R, then 
LT(/) = LT(/''), from which the second part follows. □ 

ASSUMPTION: For the rest of §5, we fix a sandpile graph G as in Definition |5.8[ and 
a sandpile monomial ordering on R. We assume the vertices are numbered so that 
Xi > Xj if i < j. 

The utility of a sandpile monomial ordering becomes apparent when one consid- 
ers topplings of sandpiles. 

Proposition 5.10. Let a, 6 G NV be distinct configurations on G such that a b, 
i.e., b is obtained from a by a sequence of vertex firings. Then, > x'^ with respect 
to the sandpile monomial ordering we have fixed on R. 

Proof. Each vertex firing deceases the size of the corresponding monomial. The 
reason is that either the vertex firing shoots sand into the sink, decreasing the total 
degree of the corresponding monomial, or it shoots sand to a vertex closer to the 
sink, in which case the corresponding monomial has more of the later indetermi- 
nates. □ 

We now proceed to compute a Grobner basis for the toppling ideal. Let 

E:ZV R 

£ ^ x'^~x'\ 

Then define T = E o A : ZV — ?► R. Thus, T{vi) is the i-th toppling polynomial, 
defined earlier, and for any configuration c, we have x'^%T(vi) — x'^ where c' is 
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the configuration obtained from c by firing Vi until Vi is stable. Morever, if cr is a 
firing-script, then x"^ % T(cr) yields the monomial corresponding to the configuration 
formed by firing a as many times as legal from c. The following theorem appears 
in the Bachelor's thesis of the second author, 27 . 

Theorem 5.11. Let b be a burning configuration, and let erf, be its script. Then 

Vt = {r{a) : < tr < (Tb} 

is a Grobner basis for I{G). 

Proof. We have im(T) C I{G) by definition of /(G). On the other hand, T{vi) is 
the j-th toppling polynomial and T{<J}y) x'' — 1. So I{G) — Span£{im(T)} by 
Proposition |4.2| 

We need to show that all S'-polynomials of Vb reduce to by F;,. Let ui and (J2 
be scripts with cti, (72 < <Jb- Write 

r((T,) = x^t - 

for i = 1,2 where c~ is the configuration obtained from by firing script ai. 
Hence, x'^'i is the leading term of T{<Ji) for each i. Define 

lcm(2:^^ , x'^i) 



X ' = 



X ' 



for i = 1, 2 so that ai + c'^ = 02 + c^ = c for some configuration c. We must show 
that the S'-polynomial, 

S{T{a,),T{a2)) = x^-Tia,) ~ x-^T{a2) 

reduces to 0. Since both scripts ai and (T2 si's legal from c, so is the script a — 
max(cri,(T2) defined by ay — max((Ti^t,, (T2,t))- Note that cr < ab. Letting c' be the 
configuration obtained by firing cr, we have the sequence of legal script- firings 



for i = 1,2, which shows that the S'-polynomial reduces to using the elements 
r(cr - CTj) for i = 1,2. □ 

Remark 5.12. In the case of an undirected graph, one may take the burning script 
to be 1, the vector whose components are all ones. Thus, the script-firings that 
are relevant in constructing the Grobner basis, described in the statement of the 
previous theorem, can be identified with firing subsets of vertices (none more than 
once). The paper [9j goes further, in this case, to describe a minimal Grobner basis, 
i.e., one in which each member has the property that none of its terms is divisible 
by the leading term of any other member. It consists of the subset of the Grobner 
basis elements described in the previous theorem corresponding to X C V such that 
the subgraphs of G induced by X and by ^ \ X are each connected. It would be 
interesting to see if this result could be generalized to the case of directed graphs. 

Theorem 5.13. Each nonnegative configuration is equivalent to a unique super- 
stable sandpile modulo L, and 

{x'^ : c is a superstable configuration) 

is the normal basis for R/I{G) with respect to the sandpile monomial ordering. 
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Proof. Two nonnegative configurations are equivalent modulo C if and only if their 
corresponding monomials are equivalent modulo the toppling ideal, I{G). In detail, 
first let ci , C2 G N" and suppose 

Ci - C2 = ^ = <'+ - ^" G £. 
Then ci > £+ and C2 > Define e = Ci - f + = ca - > 0. Then 

Conversely, suppose x'^^ —x'^^ G I{G)- Identify the sandpile group S{G) with Z"/£. 
Let 

V': C[a;i,...,x„] ^ C[Z"y£] 

be the mapping into the group algebra where is the image of the i-th standard 
basis vector for Z". Then I{G) = kcrt/j. Hence, 

In other words, Ci — C2 G C. 

Now let c be any nonnegative configuration. Since x'^ % T{a) = x'^ where c' is 
obtained by firing the script a as many times as is legal, the normal form for x'^ 
with respect to the sandpile monomial ordering is superstable. Since the normal 
form is unique, so is this superstable element. □ 

Remark 5.14. As noted in ^ we have R/I{G) « C[S{G)]. Hence, by the previous 
theorem, we see that the sandpile group can be thought of as the set of superstables 
where the sum of superstables ci and C2 is taken to be log(a;'^'^x'^^ %/(G)). 

The following can be found in Tf] for the case of Eulerian graphs. Here we 
extend the result to general sandpile graphs (for which the underlying graph is a 
directed multigraph). 

Corollary 5.15. A configuration c is superstable if and only ifc^^a-^ — c is recurrent. 

Proof. By Theorems |5.13| and |2.18[ the number of superstable configurations is 
equal to the number of recurrent configurations. Thus, is suffices to show that if c 
is superstable, then Cmax — c is recurrent. 

Let 6 be a burning configuration for G with burning script at. Since c is su- 
perstable, there exists Ui G supp(CTh) such that (c — Aai,)^ < 0. Similarly, there 
exists U2 G supp(o'f, — Ml) such that (c — A{ai, — ui))t,2 < 0. Continuing, we find a 
sequence of nonsink vertices mi, . . . , Ufc such that X)i=i ''^i ~ '^b and for 1 < j < A:, 

It follows that ui, . . . , Mfe is a legal firing sequence for Cmax— c+6, reducing Cmax—c+b 



to Cmax ~ c. Hence, Cmax — c is recurrent by Theorem 2.27 □ 



In light of Corollary |5.15[ we say that the superstables are dual to the recurrents. 
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6. Zeros of the toppling ideal 

Given any ideal / G i?. = C[a:i, . . . , x„], the set of zeros of / is 

Z{I) ^{p£ C" : f{p) = for all / G /}. 

In this section, our goal is to describe the set of zeros of the toppling ideal. 

Proposition 6.1. Let G be a sandpile graph. Then the set of zeros of its toppling 
ideal, I{G), is finite. 



Proof. Since I{G) is the lattice ideal for a square matrix of full rank, Theorem 3.1 (3| 



guarantees that the set of zeros is finite. However, we will give a direct proof. We 
have seen that 

R/I{G) « C[5(G)], 

and thus, R/I{G) is a finite-dimensional vector space over C. For each indetermi- 
nate Xi G R, consider the powers \,Xi,x1, . . . By finite-dimensionality, the image 
of these powers in the quotient ring are linearly dependent. This means there is a 
polynomial fi in one variable such that fi{xi) G I{G). Each fi will have a finite 
number of zeros, and thus, for each z, we see that the there are a finite number of 
possible i-th coordinates for any zero of the toppling ideal. □ 

Remark 6.2. In fact, the i-th coordinates of the zeros of the toppling ideal are the 
eigenvalues of the multiplication mapping 

R/I{G) R/I{G) 
9 x^g 

See [TT], for instance. 

6.1. Orbits of representations of Abelian groups. 

6.1.1. Affine case. Let {ai,...,a„} be generators (not necessarily distinct) for a 
finite Abelian group, A. Consider the exact sequence 

(6.1) O^A^Z" ^ A-^0 

Ci I ^ di 

where A is defined as the kernel of the given mapping Z" — > A. Taking duals by 
applying Homz( • ,C^) gives the sequence 

(6.2) 1 ^ A* ^ (C'')" ^A* ^ 1, 
where A* is the character group of A. 

Remark 6.3. 



(1) Exactness of (6.2) is not immediate. The exactness at A* ^ (C^)" follows 
because is a divisible Abelian group. An Abelian group B is divisible 
if for all a G i? and positive integers n there exists b G B such that nb = a. 
(For the multiplicative group , each element has an n-th root.) Applying 
Homz( • ,B) to an exact sequence of Abelian groups (Z-modules) always 
gives an exact sequence precisely when B is divisible. The proof of this, 
in general, is not immediate. However, in the case in which we are most 
concerned, the exactness is easy to estabhsh. Suppose A = S{G) is the 
sandpile group of a sandpile graph, and suppose A is the reduced Laplacian 
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lattice, £ = im(A) ^ Z". We would like to show that the natural map, 
given by composition, 

Hom (Z" , C ) ^ Hom (£ , C ) 

is surjective. Let be given. Since the reduced Laplacian has 

full rank, given v G Z", there exist unique rational numbers a£ such that 
V = ^£ ai£, with the sum going over a basis for C (say, over the columns 
of the reduced Laplacian). Then define 0: Z" ^- by }{v) = J^i <^(^)"'- 
(2) To be explicit, denote the mapping Z" A hy (j). Then part of sequence 



(6.21 is 



A* Hom(Z",C^) « (C^)" 

X ^ X°<P (x(ai),---,x(a„)). 



We get an n-dimensional representation of A*: 

p: A* ^ (C^)" ^ GL(C") 

given by 

p{x) = diag(x(ai), • ■ • , xM)- 

In other words, the choice of generators for A induces a homomorphism of A* into 
the group of invertible nx n matrices over C. (Every n-dimensional representation 
of A* over C is a direct sum of characters of A* , i.e., of elements of A** « A. So 
this section can be regarded as saying something about representations of A* , in 
general.) 

For each z G C", define the orbit of z under p to be 

Op{z) = {p{x)z : X e = {{x{ai)zi. ■ ■ ■ .X{an)zn) ■ X € A*}. 

We will assume that no coordinate of z is zero, in which case by scaling coor- 
dinates of C", we may assume for our purposes that z = (1, . . . , 1). Thus, we are 
interested in the orbit of the all-Is vector: 

O = {p{x) ■■ X e A*} - {(x(ai), . . . , xM) : X e A*}. 

Definition 6.4. Let / C i? be an ideal. Le i?<rf denote the vector space of poly- 
nomials in R of degree at most d, and let I<d be the subspace / n R<d- The ajfine 
Hilbert function of R/I is 77 : N — N, given by 

H{d) := dime R<d/I<d = dime R<d - dime I<d- 

Theorem 6.5. Let R = C[xi, . . . , a;„] and consider 

I = {feR: f(0) = 0}, 

the ideal of polynomials vanishing on the orbit. Then 

(1) I = I{A) = Spane{a:" - : u ^ v mod A}; 

(2) The affine Hilbert function of R/I is given by 

H{d) — ^ "-i^i : n-i > for all i and Ui < d} . 

Proof. This proof is due to the first author and Donna Glassbrenner. It appears 
in [7]. Consider the matrix M^'^^ with rows indexed by A* and columns indexed by 
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the monomials of R<d (arranged in lexicographical order so that M^''-' is naturally 
nested in Af('*+^)) given by 



4=1 



Using the isomorphism 



A ^ A* 



where a(x) '■— x(fl)i can write 

n 

^x'i"=n«r(x)=a"(x) 
1=1 

where a" := n"=i € ^**. The a;"-th column of M^'') has entries a"(x) as x 
varies over A* . In other words, it is the list of all values of the function a". Thus, 
at least as far as linear algebra is concerned, the a;"-th column is a". Since distinct 
characters are linearly independent, it follows that any linear dependence relations 
are the result of columns that are equal. 

Now, the a;"-th and x"-th columns of M^''^ are equal exactly when a" = a" are 
equal. This occurs exactly when u^Oi — WiCti, which we write as (u — w) - a = 
where a :— (ai, . . . , a„). In light of exact sequence (6.1 1, this condition is equivalent 
to u — u e A. 

A vector (a„) € kerM^''' if and only if 



u i—1 



for all X G Thus, € kerM'^''^ if and only if the polynomial p = a„x" 
vanishes on O, i.e., p G I. Thus, elements of I<d correspond exactly with linear 
combinations among the columns of M^''^ . As these relations are due to equality 
among columns, as already noted, part [l] follows. For part [2j note that we have 
just shown that 

dim/<rf = dimi?<d - rank M^'*^ 
Since distinct characters are linearly independent, 

rankM*-''-' = # {J2^=i ''^i^i ■ n-i > for all i and J2i ^ d} . 

□ 

Returning to the case of the toppling ideal, the exact sequence 
^ Z" A Z" ^ S{G) 



has the form of exact sequence (6.1). The generators are the configurations 
having exactly one grain of sand. 
Corollary 6.6. 

(1) The toppling ideal is the set of polynomials vanishing on an orbit O of a 
faithful representation ofS{G)*. 

(2) The set of zeros of the toppling ideal is the finite set, O. It thus has the 
symmetry of S{G)* , which is isomorphic to the sandpile group. 
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(3) // Hq is the affine Hilbert function for the toppling ideal, then Hc^d) is 
the number of elements of '27^ / C represented by configurations containing 
at most d grains of sand. It is thus the number of superstable configurations 
of degree at most d or, eguivalently, the number of recurrent configurations c 
such that 

deg(c) > deg(cmax) - d. 

Proof. Part ([T]) follows directly from the first part of Theorem |6.5[ For part ^ , 
since O is a finite collection of points in C", and I{G) = I{0), it is a basic result 
of algebraic geometry that the set of zeros of /(G) is O. Part (|3| is immediate 
from the second part of Theorem |6.5| and the fact that r is recurrent if and only if 
Cmax — J' is superstable. □ 



Remark 6.7. Note that part ([s]) also follows directly from Theorem 5.13 



6.1.2. Projective case. An ideal J in S = C[xi, . . . , a;„_|_i] is homogeneous if it has a 
set of homogeneous generators. The set of zeros of J is a subset of projective space: 

Z{J) = {peF" : f{p)=0 for all homogeneous / G J}. 

The ring S/J is graded by the integers: {S/J)d ■= Sd/Jd- 

Definition 6.8. The Hilbert function of S/J is i7 : N ^ N, given by 

H{d) := dimciS/J)d. 

Continuing with the notation from |6.1.1[ define the homogenization of A as 

Consider the exact sequence 



^ A'' ^ Z"+i A A e Z ^ 0, 



where 



M 



fli ... a„ 
1 ... 1 1 



Apply Hom( • , C^) to get 

l-^A*xC'^ ^ (C'^)"+i ^ (A'')* ^ 
{X,z) ^ (x{a'i)z,...,x{a'n)z,z) 
and the corresponding representation 

A* y-C"" GL(C"+i) 

(X,z) ^ diag(x(ai)z,...,x(a„)z,z). 
The orbit of (1, . . . , 1) under this representation is 

= {(X(ai), . . . , X(a„), 1) e P" : X e C P". 
Thus, is the projective closure of the orbit O from the previous section. 

Theorem 6.9. Let — (ai, . . . , a„, 0). 

(1) The homogeneous ideal defining is the lattice ideal for A''^ , the homoge- 
nization of the lattice ideal for A; 

= {a;" - : u = w mod A''}. 
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(2) The Hilbert function for (i.e., the Hilbert function of S/l'^) is 

H{d) ^#{h-a^ eA:be N"+^ and deg(6) = d}, 

which is the same as the affine Hilbert function for O. 

Proof. Since O'' is the projective closure of O, its ideal is l'\ the homogenization 
of the ideal defining O, which is given by {cc" — x"" : u = v mod A''}. The second 
part of the theorem follows from part [2] of Theorem 6.5 and the isomorphism of 
vector spaces 

Sd — > R<d 
f ^ /U„+i=i, 

with inverse g{xi, . . . ,a;„) i-^ x^_^-^ g{xi/xn+i, ■ • . ,x„/a;„+i). □ 

Corollary 6.10. Suppose A — C, the reduced Laplacian lattice of G, and that 
/S.(vn+i) G Span2{A(tii) : 1 < i < n} so that = C, the full Laplacian lattice (see 



the comments preceding Proposition 4-- 8) 



(1) The homogenization of the toppling ideal is the ideal generated by all homo- 
geneous polynomials vanishing on an orbit of a faithful representation 
of{Z"+^/C)*. 

(2) The set of zeros of the homogenization of the toppling ideal is the finite set 

having the symmetry of S{G)* . 

6.1.3. The h-vector. Let AHc; denote the first differences of the affine Hilbert func- 



tion of a sandpile graph G. So AHG{d) :— Haid) — Hcid — 1). By Theorem 5.13 
the value of AHcid) is the number of superstable configurations of degree d. 

Definition 6.11. Let hd AHoid). The postulation number for G is the largest 
integer i such that hi ^ 0. The ft,- vector for G is /i = (/iq, ■ ■ • , hi). The Hilbert- 
Poincare series for G is Pciv) — J2i=o ^iV^- 

Example 6.12. Continuing Example |2.15[ the /i-vector for the sandpile graph in 
Figure [l] is (1,3,6,7,4). 

Let the vertices of G be {vi, . . . ,Vn+i\ with as the sink, as usual. Let 

l'^ ^ Ih C S = C[a;i, . . . , Xn+i] be the homogenization of the the toppling ideal and 
the homogeneous toppling ideal for G, respectively. These two ideals are identical 
when the hypothesis of Proposition |4.8| is satisfied. In any case, their zero-sets 
satisfy Z{I^) D Z{Ih). Pick a linear polynomial f € S that does not vanish at any 
point of Z[I^). For instance, we could take f ~ Xi for any i. Multiplication by / 
gives rise to the commutative diagram with exact rows 

> S/I^ — ^ S/I^ > S/{I^ + (/)) > 



> S/h s/h > s/ih + if)) > 0. 



By this diagram and Theorem |6.9[ we have the following relations among the first 
differences of Hilbert functions: 

(6.3) hd = AHaid) = AH s/ 1, id) = Hsf^i.+^f)){d) 

>HsninHf))id)^^Hs,iM- 
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6.1.4. The Tutte polynomial. Now let G — {V, E) be any (weighted, directed) graph, 
and e e -E. Let G — e denote the graph obtained from G by replacing wt(e) by 
wt(e) — 1. In other words, imagine the endpoints e~ and e"*" attached by wt(e) 
edges, and remove one of these edges to obtain G — e. In particular, if wt(e) = 1, 
this amounts to removing the edge e. Let G/e denote the graph obtained from G 
by identifying the endpoints of e and lowering the weight of e by one. We refer to 
these two operations on G as deletion and contraction. The edge e is called a bridge 
if G — e has more components than G. 

Definition 6.13. Let G be an undirected, weighted graph. Define the Tutte poly- 
nomial, TQ{x,y) for G recursively, as follows. If E consists of i bridges, j loops, 
and no other edges, then 

TG{x,y) := x'y^. 

In particular, = 1 if G has no edges. Otherwise, if e G E' is neither a bridge nor 
a loop, then 

Tq ■= Tc-e + Ta/e- 

It turns out the the Tutte polynomial is well-defined, independent of choices for 
deletions and contractions. It is well-known that 

Gg(x) = (-1)#^-«(«)x''«=)Tg(1 - X, 0), 

where Cg is the chromatic polynomial of G and k{G) is the number of components 
of G. The following result relates other specializations of the Tutte polynomial to 
the algebraic geometry of sandpiles. 

Theorem 6.14 (Merino [21] )• Let G be an undirected sandpile graph with postu- 
lation number i. Then 



TcH^y) =Y.h,_.,y\ 

Corollary 6.15. Let G be as in Theorem \6.14\ Then 

(1) the Hilbert-Poincare series for G is y^ T{1, 1/y); 

(2) if d is the degree of the maximal stable configuration on G, then y"^^^ T{l,y) 
is the generating function for the recurrent configurations of G (by degree); 

(3) Tq{1, 1) is the order of the sandpile group of G; 

(4) Tg'(1,0) is the number of maximal superstable (or the number of minimal 
recurrent) configurations ofG. 



Proof. These results follow immediately from Theorem 6.14 Part ^ uses the fact 



that c is superstable if and only if Cmax — c is recurrent. □ 

Example 6.16. Figure [7] shows the construction of the Tutte polynomial of a 
graph G. We have T(l, y) = 4 + 3y + and T{1, 1) = 8. Fixing the southern-most 
vertex of G as the sink gives a sandpile graph with /i-vector (1,3,4) and sandpile 
group of order 8. 

6.1.5. Cayley-Bacharach property. Let X C P" be a finite set of points in projective 
space, and let I{X) <z S :— C[xi, . . . , Xn+i] be the ideal generated by the homoge- 
neous polynomials vanishing on X. If Hx is the Hilbert function of S/I{X), then 
Hx{d) is the number of linear conditions placed on the coefficients of a general 
homogeneous polynomial of degree d in 5 by requiring the polynomial to vanish on 
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Taix, y) ^ X + 2x'^ + + [1 + 2x)y + y^ 



Figure 7. The Tutte polynomial of G. 



the points of X. Thus, Hx is a monotonically increasing function which is eventu- 
ally constant at \X\. The first value at which Hx takes the value |^| is called the 
postulation number for X. 

Definition 6.17. A finite set of points X C P" is Cayley-Bacharach if it satisfies 
one of the following equivalent conditions. 

(1) For each p G X, and for each d <E N, 

Hx\{p}id)=mm{Hxid),\X\-l}. 

(2) Every homogeneous polynomial with degree less than the postulation num- 
ber for X and vanishing on all but one point of X must vanish on all of X. 

Proposition 6.18. The set of zeros of the homogeneous toppling ideal is Cayley- 
Bacharach. 

Proof. By Proposition 1.14 of [IS], for any finite set of points, X, there is always at 
least one point p for which condition ([T]) of Definition 6.17 holds. However, in our 
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case, X is the orbit of a linear representation of the sandpile group. Thus, given 
any two points p,q £ X, there is a hnear change of coordinates of P" sending p 
to q. A hnear change of coordinates does not change the Hilbert function. Hence, 
condition ([T]) holds for all points oi X. □ 

Remark 6.19. Let X be the set of zeros of a homogeneous toppling ideal and define 
the first differences of its Hilbert function by AHx{d) = Hx{d) — Hx{d — 1) for all 
d G Z. It follows from results in [1^] and the fact that X is Cayley-Bacharach, that 
if the last nonzero value of Ax is to, then there is a collection of m points Y <Z X 
such that X \ y is Cayley-Bacharach. Moreover, if m = 1, then every subset of X 
of size \X\ — 1 is Cayley-Bacharach. 

7. Resolutions 

In this section, we consider the minimal free resolution of the homogeneous top- 
pling ideal, summarizing some of the results in |32) . For further work on resolutions 
of toppling ideals, see [20]. First, we recall the language of divisors on graphs from 
[3] (extended to directed multigraphs) . Let G be a directed multigraph as in S|2] 
The free Abelian group "LV on the vertices of G is denoted div(G'), and its el- 
ements are called divisors. The degree of a divisor D = ^y^y Dy v G div(G'), is 
deg(_D) := J2vev ^ divisor is principal if it is in the Laplacian lattice £, defined 
in fj2] Divisors D and D' are linearly equivalent, written D ^ D' , ii D — D' is prin- 
cipal. Note that linearly equivalent divisors must have the same degree. The group 
of divisors modulo linear equivalence is the class group of G, denoted C1(G). In the 
case where G is an Eulerian sandpile graph, using the notation of Proposition |2 . 3 1 [ 
there is an isomorphism 

C1(G) ^ Z ® ZVa/C 

D ^ {deg D,D~ {deg D)s), 

where ZVq/C is isomorphic to the sandpile group S{G). 

We will usually denote a divisor class [D] e C1(G) by just D, choosing a repre- 
sentative divisor, when the context is clear. A divisor D — '^y^y Dy v is effective 
ii D > 0. The collection of all effective divisors linearly equivalent to D is called 
the (complete) linear system for D and denoted \D\; it only depends on the divisor 
class of D. The support of a divisor D is supp(£') :^ {v E V : Dy ^ 0}. 

One might think of a divisor as an assignment of money to each vertex, with 
negative numbers denoting debt. Just as with configurations in the sandpile model, 
the Laplacian determines firing rules by which vertices can lend to or borrow from 
neighbors. Two divisors are linearly equivalent if one can be obtained from the 
other by a sequence of vertex lendings and borrowings. The complete linear system 
corresponding to a divisor is nonempty if there is a way for vertices to lend and 
borrow, resulting in no vertex being in debt. 

7.1. Riemann-Roch. To recall the graph-theoretic Riemann-Roch theorem of [3], 
let G = {V, E) be an undirected graph. Define the genus of G to be 

Define the dimension of the linear system \D\ for a divisor D on G to be 
r{D) max{A: £ Z : - 7^ for all > with deg(£') = fc}. 
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with r{D) := —1 if \D\ = 0. Note that r{D) depends only on the divisor class of D. 
Define the maximal stable divisor, 



-Dniax ^(deg(-y) - l)v, 



and the canonical divisor, 



Theorem 7.1 (Riemann-Roch Theorem j3j). Let G he an undirected graph. For 
all D e div(G), 

r{D) - r{K - D) = Aeg{D) + l-g. 

Remark 7.2. This Riemann-Roch theorem is generalized in [5D] to the context of 
certain monomial ideals, relating it to Alexander duality in combinatorial commu- 
tative algebra. From that point of view, the relevant monomial ideal for us is the 
ideal generated by the leading terms of a homogeneous toppling ideal with respect 
to a sandpile monomial ordering. It is noted that these monomial ideals are studied 
by Postnikov and Shapiro in .28j . 

7.2. Resolutions and Betti numbers. Let G be an arbitrary directed multi- 
graph. Identify the vertices of G with the set {1, . . . , n -I- 1}, with n + 1 being 
the sink. The polynomial ring S — C[xi, . . . , a;„+i] is graded by the class group 
by letting the degree of a monomial be D G C1(G). For each D e C1(G), let 
Sd be the C-vector space generated by the monomials of degree D, and define the 
twist, S{D), by letting S{D)f := S'(d+f) for each F e C1(G). 

Let / := Ih{G) be the homogeneous toppling ideal. A free resolution of / is an 
exact sequence 

^ / <^ Fi 4^ F2 ^ • • • 4^ ^ 0, 
where each Fi is a free Cl(G)-graded 5- module, i.e., 

Deci(G) 

for some nonnegative integers /3i,_D, and where each (j> preserves degrees. The length 
of the resolution is r . A free resolution is minimal if each of the Pi^D is the minimum 
possible from among all free resolutions of /. In this case, the fii^D are called the 
Betti numbers of /. For instance, Pi^d is the number of polynomials of degree D in 
a minimal generating set for /. We also define the z-th coarsely graded Betti number 
of / by ft = Y.D&C\{D)P%D- 

The following theorem states a well-known fact about resolutions of sets of points 
in projective space (the Cohen-Macaulay property). 



Proposition 7.3. The length of the minimal free resolution of the homogeneous 
toppling ideal is n, the number of nonsink vertices. 
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Vl 




V2 



V3 



1 

V4 



Figure 8. A Gorenstein sandpile graph G with sink V4. 



Example 7.4. Let G be as in Figure [s] and let / = 1(G)''. Then 



0^ 



0101 
0110 
1010 
1100 
0011 



1011 
1101 
1110 
2010 
0111 



is a minimal free resolution for /, where the are given by 



X1X3 X1X2 



- [ 




X2X4 X2X3 - 


XIX4 


2 

X2 ~ 






Xi 





X4 


" 






-X2 


Xi 





—X4 


1 = 


X4 


X3 





Xi 













-X3 


-X2 


Xi 










X4 


X3 


-X2_ 




'- xl- 


- X3X4 










-X1X2 + x\ 








2 = 


X2 
X2X3 

X3 


+ X1X3 

— X1X4 

- X2X4 









X3X4\ 



The grading of the S'-modules is indicated below each of them. For example, the 
last ^-module is ^(-(l, 0, 2, 2)). 

The Betti numbers of / may be understood topologically. For D G C1(G), 
define the simplicial complex A^i on the vertices of G by G Ajj if and only 
ii W C supp(i?) for some E € \D\. The following version of Hochster's formula 
appeared as Lemma 2.1 of |25| . 

Theorem 7.5. The Betti number /3i.£> is the dimension of the [i — l)-th reduced 
homology group Hi^i{l^D\<C) as a C-vector space. 



Example 7.6. Let G again be as in Figure |8] For D = vi 
Example 7.4 that (32, d = 1- We have 



"^3 + ^4, we saw in 



\D\ = {D,V2 + 2v3,3vi,2v2 + V4}, 

so the simplicial complex A/j is as pictured in Figure [g] Note dime i?i ( Ad ;C) = 1, 
as asserted by Hochster's formula. 
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Vl 




Figure 9. The simplicial complex Ad for Example 7.6 



7.3. Minimal recurrents. Again specialize to the case of an undirected graph G. 
As part of the Rieniann-Roch theory, one defines the non-special divisors on G to 
be 

JV — {De div(G) : deg(i:>) =g-l and \D\ = 0}. 

By the Riemann-Roch theorem, if deg{D) > g—1, then |D| 7^ 0. So the nonspecial 
divisors are the divisors of maximal degree having empty linear system. 

Fix s G V and consider the sandpile graph G — {V,E,s). A recurrent configu- 
ration c on G is minimal ii c — v is not recurrent for any nonsink vertex v. It is 
well-known that (since G is undirected) the minimal recurrent configurations are ex- 
actly the recurrent configurations of minimal degree, namely of degree #£' — deg(s). 



(This result follows from Dhar's algorithm (cf. j |2.4| and the proof of Theorem 8.28[ ).) 
Similarly, one says that a superstable configuration c is a maximal if c + v is not 
superstable for any nonsink vertex v. By Corollary |5.15[ it follows that the maximal 
superstable configurations are exactly those of degree g. 

We say that a divisor Z? on G is unstable if D^, > deg(u) for some v G V and 
that D is alive if there is no stable divisor in \D\. Further, D is minimally alive if for 
all w € 1^, we have that D — v is not alive. It is shown in f32| that a divisor D is alive 
if and only if D ^ c + k s for some recurrent configuration c and some k > deg(s), 
and D is minimally alive if and only if 13 ^ c + deg(s) s for some minimal recurrent 
configuration c. 

It is shown in 3 that a set of representatives for the distinct divisor classes of 
the non-special divisors is 

{c — s : c a maximal superstable configuration}. 

Thus, the non-special divisor classes are given, essentially, by the maximal G- 
parking functions. 

Suppose that v is a nonspecial divisor. We may assume v = c — s for some 
maximal superstable configuration c. Then 

-Dmax (c,„ax - c) + deg(s) S. 

Since Cmax — c is a minimal recurrent configuration, i'max — Jy is minimally alive. 
Similarly, one may show that if D' is a minimally alive divisor, then -Dmax — D' is 
nonspecial. Thus, on an undirected graph there is a bijective correspondence be- 
tween: minimal recurrent configurations, maximal superstable configurations, max- 
imal G-parking functions, acyclic orientations with s as the unique source vertex 



(cf. Theorem 2.25), minimally alive divisors, and non-special divisors. In particular, 
the cardinality of these sets does not depend on the choice of sink. 

The following is Theorem 3.10 of [35]. The proof is included here for the sake of 
completeness. 
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Theorem 7.7. Let G be an undirected graph and D E div(G). Let r = — 1, 

the length of a minimal free resolution for G. Then the highest nonzero Betti 
number, 13^, is the number of minimal recurrent configurations on G. We have 

f3rM + 

if and only if D is minimally alive (in which case deg{D) = jj=E). 

Proof. First note that by Theorem 1 7. 1[ a divisor v is nonspecial if and only \i K --v 
is nonspecial. Indeed, if deg(i/) = g — 1, then 

deg(if -v)^{2g~2)~g~-l^g-l 



so that Theorem 7.1 gives r{K — v) = r{v). 

By Theorem 7.5 we have j3r,D — dime i?r-i(AD; C). Since for any D E div(G') 
the simplicial complex Ajj has vertices, ^ if and only if Ajj is the 

boundary of an r-simplex. Thus, ^ if and only if f3r,D = 1, or equivalently: 
(i) no G 1 13 1 has full support, and (ii) for every v E V there is some E E \D\ with 
C supp(£;). 

Suppose D is minimally alive. Then 13max — 13 is nonspecial by the discussion 
preceding the statement of the theorem. Let — K — (i^max ~ D), so that v is 
also nonspecial. In particular, |_D — 1| = \v\ = 0, so no divisor E E \D\ has full 
support. Now Hx V E V and let F = D — 1 + v. Note that a divisor EE \D\ 
satisfies V \ {v} C supp(£') if and only if £■ — 1 + u E \F\. So to complete the 
proof that /J^.d 7^ it suffices to show that |F| 7^ 0. Note that deg(F) — g. 
Since K — F + v = D^jii^^ — D, we have K — F + v nonspecial, and it follows that 



r{K — F) = —1. Thus, by Theorem 7.1 we have r{F) = as desired. Hence, D 
satisfies (i) and (ii). 

On the other hand, suppose D satisfies (i) and (ii) above, and let i/ = D — 1. 
Then = follows from (i), and therefore -Dmax — is alive. On the other hand, 
for every v E V we have from (ii) that |^ + u| 0, whence (-Dmax — ly) ~ v is 
not alive. Thus, -Dmax — v is minimally alive, so that i' is nonspecial. But then 
K — v = -Dmax — -D is also nonspecial, implying D is minimally alive. □ 

Example 7.8. We summarize many of the results of this paper using the graph G 



of genus g — 2 in Figure 10 The mathematical software Sage |30j was used for 



some of the calculations. The sandpile group for G is cyclic of order 8. Its toppling 




Figure 10. Genus two graph G. 

ideal is / = [x^ — yz, y^ — xz, — xy, yz — 1), and its homogeneous toppling ideal 
is 

T rh f 2 3 3 22222\ 

Ifi — ^ — ~ yz, y — xzs, z — xys, yz — s , xz — y s,xy — z s). 
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Letting w = exp(27ri/8), the zeros set of / is 

Z{I) = {{{-ly : < J < 7} C C3, 

which forms a cychc group of order 8 under component-wise muhipUcation. With 
respect to the sandpile monomial ordering (grevlex) for which x > y > z > the 
normal basis for the coordinate ring of G is the spanned by 8 monomials: 

R/I = <C[x, y, z]/I = Span{l, x, y, z, xy, xz, 2/^ z^}. 

The exponent vectors of the normal basis give the superstable configurations: 

(0, 0, 0), (1, 0, 0), (0, 1, 0), (0, 0, 1), (1, 1, 0), (1, 0, 1), (0, 2, 0), (0, 0, 2), 
and dualizing, c — )■ Cmax ~ c, gives the recurrent configurations: 

(1, 2, 2), (0, 2, 2), (1, 1, 2), (1, 2, 1), (0, 1, 2), (0, 2, 1), (1, 0, 2), (1, 2, 0). 

(We use the notation (ci, C2, C3) := ci x + C2 y + C3 z.) 

From the degrees of the monomials in the normal basis, one sees that the affine 
Hilbert function for G is 

i7G(0) = l, Hg(1) = 3, -ffG(2)=4 

with postulation number 2 (equal to 5, the degree of the maximal superstables). 
The Tutte polynomial for G was calculated in Figure [7j and in accordance with 
Corollary |6.15[ the Hilbert series for G is 

y'^TailMy) - l + 3y + V. 

The minimal free resolution for G is 

00 



0< /f- 





0. 


SU 


Olio 


0121 


0122 


2000 


0211 


0212 


0030 


1201x2 


1022 


1020 


1120 


1202 


0300 


1021x2 




1200 


1210 






0220 





The Cl(G)-degrees are listed in x, y, z, s order. The degrees of the highest nonzero 
Betti numbers correspond to the minimal recurrent configurations as prescribed by 
Theorem |7.7| For instance, the degree 0122 corresponds to the minimal alive divisor 
y + 2z + 2s and to the minimal recurrent configuration (0, 1, 2). Thus, ~ Ha{2), 
and the degrees of each of these divisors is 5 = ^E. 

As an example of Hochster's formula (Theorem 7.5 1, let D — 1021 — x + 2z + s. 
The complete linear system for D is 

l^l = {1021,2200,0202,0310}, 

and Ajj is the simplical complex pictured in Figure |11[ We have 

/32,i02i =dimci/i(1021;C) = 2. 
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Figure 11. The simplicial complex for Example 1 7. 8 



7.4. Conjecture. Let G — {V, E, s) be an undirected sandpile graph. For U CV, 
let G\u denote the subgraph of G induced by U, i.e., the graph with vertices U 
and edges e £ E such that both endpoints of e are in U. A connected k-partition 
or k-bond of G is a partition 11 = U^Lj^Vi of V such that Gly. is connected for all i. 
The corresponding k-partition graph, Gu, is the graph with vertices {^i, . . . ,14} 
and with edge weights 

wt(Vi, Vj) — #{e e E : one endpoint of e is in Vi and the other is in Vj}. 

We consider Gu to be a sandpile graph with sink vertex Vi, where i is chosen so 
that s &Vi. 

The following conjecture appears as Corollary 3.29 in [32] ■ Using the mathemat- 
ical software Sage, it has been verified for all undirected, unweighted graphs with 
fewer than 7 vertices. 

Conjecture 7.9. LetVk denote the set of connected k-partitions ofG. Then 
Pk — #{c : c a minimal recurrent configuration on Gn}- 



Example 7.10. Figure 12 displays the 5 connected 3-partitions of G along with 
their corresponding 3-partition graphs and /i-vectors. The top value of each h- 
vector is the number of minimal recurrent configurations (or maximal superstable 
configurations) on the partition graph. Summing these top values gives /32 for G. 



As a corollary to Conjecture 7.9 it is shown in [32 that 



Corollary 7.11. // Conjecture \ 7.SI\ is true, then the following five statements are 
also true. 

(1) The number of polynomials, Pi, in a minimal generating set for the homo- 
geneous toppling ideal of G is equal to the number o/cuts (i.e., the number 
of connected 2-partitions) of G. 

(2) For a tree on n vertices, /3fc = ("^ ) . 

(3) // the weight of an edge of G is changed from one nonzero value to another, 
the j3k do not change. 

(4) // G' is obtained from G by adding an edge to G (between two vertices of 
G), then f3k{G) < /3fc(G') for all k. 

(5) For the complete graph on n vertices, K^, we have that 13^. is the number 
of strictly ascending chains of length k of nonempty subsets of [n — 1] :— 
{l,...,n~l}. 



Remark 7.12. Corollary 7.11| ([sj) is proved, independently, in 
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connected 
3-partitions 




i 



i 



i 



i 



i 



T 



3-partition 
graphs 




/i-vector 



1 2 2 



1 2 2 



1 2 1 



1 2 2 



1 2 2 



Figure 12. Second Betti number: /32 = 2 + 2 + 1 + 2 + 2 = 9. 



8. GORENSTEIN TOPPLING IDEALS 



This section characterizes toppUng ideals that are complete intersection ideals 
and gives a method for constructing Gorenstein toppling ideals. 

8.1. Complete intersections. If y C P" is the solution set to a system of homo- 
geneous polynomials, then y is a complete intersection if the ideal generated by all 
homogeneous polynomials vanishing on V can be generated by a set of polynomials 
with cardinality equal to the codimension of V in P". Specializing to the case of 
sandpiles, we get the following definition. 

Definition 8.1. Let G = (V, E, s) be a sandpile graph with homogeneous toppling 
ideal /. Then G is a complete intersection sandpile graph if / is generated by \V\ — 1 
homogeneous polynomials. (We also say that / or the set of zeros of / is a complete 
intersection.) 

Let £ be a submodule of Z"+^ of rank n whose lattice ideal /(£) is homoge- 
neous. Then /(£) is a complete intersection if it is generated by n homogeneous 
polynomials. 

Remark 8.2. The lattice ideal I{C) is generated by homogeneous polynomials if 
and only if deg(w) := J2i = ^ for all w £ £. 

For the following, recall from S|2]that sandpile has an absolute sink if its sink has 
outdegree 0. 

Definition 8.3. For i — 1,2, let Gi = {Vi,Ei,Si) be a sandpile graph with edge- 
weight function wt^ and absolute sink Sj. Suppose that the two graphs are vertex- 
disjoint. Let G be any graph with vertex set V — Vi U V2, and edge-weight func- 
tion, wt, satisfying the following 

(1) wt(e) = wti(e) if e G Ei, 

(2) wt(e) = wt2(e) if e e E2, 

(3) wt(u, w) = if (u, v) g (Vi X V2) U {V2 X Vi), 

(4) wt(si,?;) > for some v £ ¥2- 
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We consider G to be a sandpile graph with S2 as its absolute sink. Let A := Aq 
be the Laplacian of G, and define 

D A(si)|y, = ^ A{si%v, 
veVi 

a divisor on Gi. Then G is a wiring of Gi into G2 with wiring divisor D if \D\ ^ 0, 
i.e., if the complete linear system for 13 as a divisor on Gi is nonempty (cf. Sj?]). 

Thus, to form a wiring of Gi into G2, one connects si into G2 with at least 
one edge and then adds edges from si back into Gi as determined by a divisor, D, 
on Gi having a nonempty complete linear system. There always exists some wiring 
of Gi into G2. For instance, we could take D — k si for any fc > by connecting Gi 
to G2 with k edges from si into G2 (and no edges from si back into Gi). 




Figure 13. A wiring of Gi into G2. 



Notation 8.4. For any sandpile graph G = {V,E,s), with Laplacian Ag, we let 
Aq ~ AgI^. Thus, A^ : ZV — >■ ZV, and in terms of matrices, Aq is formed 
from Ag by removing the column corresponding to the sink — a column of zeros 
if G has an absolute sink. We will call A^ the restricted Laplacian of G. 

With this notation, if G is a wiring of Gi into G2, then 

A^ a \ 

where exactly one entry of a is positive (corresponding to si) and /3 < 0. The last 
column corresponds to si, and the wiring divisor is D — a. 

If Gi is a single point with no edges, then we regard A^ as the 1x0 empty 
matrix, and a will be a single integer, as in the following example. 




Example 8.5. Let Gi be the graph with a single vertex si and no edges. Let G2 
have vertex set {w2,W3,S2} and edge set {(W27S2): {^3,32)}- Figure 14 illustrates a 
wiring, G, of Gi into G2. The wiring divisor is 13 = 2si. The restricted Laplacian 
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V2 




V3 



Figure 14. The wiring G for Example [83 
of G is, with respect to the indicated vertex ordering, 





V2 


V3 






( ' 







V2 


1 





-1 


V3 





1 


-1 


S2 




-1 





Definition 8.6. A directed multigraph G is completely wired if it is a single vertex 
with no edges or if it is the wiring of one completely wired graph into another. 

Example 8.7. Every directed acyclic graph is completely wired. 

Definition 8.8. An integral matrix is mixed if each column contains both positive 
and negative entries. An integral matrix is mixed dominating if it does not contain 
a mixed square submatrix. 

Empty d X matrices are mixed dominating by convention. The following two 
theorems are established in ^3] and [T3]. 

Theorem 8.9. Let C be a submodule of Z"+^ of rank n such that the associated 
lattice ideal I{C) is homogeneous. Then I{C) is a complete intersection if and only 
if there exists a basis ui, . . . , w„ for C such that the matrix whose columns are the Ui 
is mixed dominating. 

Theorem 8.10. If M is a mixed dominating matrix, then by reordering its columns 
and rows we may obtain 



M' = 



where the Mi are mixed dominating, a > 0, and (3 < 0. 



( Ml 







\ 


M2 





It is allowable for the matrix Mi in Theorem |8.10| to be the empty d x matrix, 
in which case we would have 



M' 





M2 



where the upper-left block is a zero matrix with d rows. A similar statement holds 
if M2 is the d X matrix, in which case we would have a lower-left zero matrix 
block. 

We now characterize complete intersection sandpile graphs. 
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Theorem 8.11. Let C be a suhmodule ofE"^^ of rank n such that the associated 
lattice ideal I{C) is a complete intersection. Then there exists a completely wired 
graph G whose Laplacian lattice is C, and hence, I{C) = I{G)^ , where I{G)^ is the 
homogeneous toppling ideal of G. 

Proof. We proceed by induction, the case n — being trivial. Let ui, . . . , m„ be a 
basis for £, and let M be the matrix whose columns are the Ui. By Remark |8.2[ 
we have deg(u.i) = for all i. (Here, deg{ui) denotes the degree of Ui as a divisor, 
i.e., the sum of the components of Ui.) By Theorems 8.9 and 8.10 we may assume 
that 

■ Ml a 
Ah P 

where the Mi are mixed dominating, a > 0, and /3 < 0. Each column of Mi and M2 
has entries that sum to zero. By our rank assumption, it follows that Mi and M2 
are matrices of full rank, each with one more row than column. By induction, 
there exist completely wired graphs Gi and G2 such that im(AG. ) — im(Mi) for 
i = 1,2. Let si be the sink of Gi. Let c be any nonnegative configuration on Gi 
with full support and contained in im(A(3j), the reduced Laplacian lattice for Gi. 
For instance, we could take c = 5 — 5° where 5 = '}2iv(zv^{ovLt(ieg{v) + 1) v. Define 
the divisor D ~ c — deg(c) si e im(AGi) = m\{/^°Q_^). Take fc e N such that 
fc • c + deg(a) si > a. Now 



M = 



M' 











A? 



a-kD 



has the same column span as M, and M' — Aq where G is the wiring of Gi into G2 
with wiring divisor a — kD. Then G is completely wired and, up to an ordering of 
its vertices, its full Laplacian lattice is C. □ 



Example 8.12. The graph of Example 4.9 is a complete intersection sandpile 
graph. It is not completely wired, but its Laplacian lattice is the same as that for 
the completely wired graph consisting of a single directed edge connecting vi to V2. 

Theorem 8.13. // the graph G is completely wired, then I{G)^ is a complete 
intersection. 

Proof. If G has only one vertex, then /(G) = {0} is a complete intersection, so we 
will again proceed by induction, now on |F(G)|. Assume |1^(G)| > 1 and that G 
is the wiring of some graph Gi with sink s into another graph G2 with wiring 
divisor D. Let /3 = Ag{s)\v2- Then 



A° — 







A5 



D 



By Theorem 8.9 and induction, there exist Mi and M2 with im(Mi) 
i = 1,2, and E e \D\, such that 



im(AG.) for 



M = 



Ml 






M2 



has the same column span as A^ and is mixed dominating. So I{Gy^ is a complete 
intersection by Theorem |8.9| □ 
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8.2. Gorenstein sandpile graphs. Having characterized complete intersection 
sandpile graphs, we proceed to give a method for constructing sandpile graphs with 
Gorenstein toppling ideals. Our basic reference for Gorenstein ideals is |15) . 

Notation 8.14. Let S = C[xi, . . . , Xn+i], and let / be a homogeneous ideal in S. 
Let Sd be the C-vector space generated by all homogeneous polynomials of degree d, 
and let /<j := / n Sd- Define A = S/I, and let Ad := iS/I)d Sd/Id- Let 

m = (xi, . . .,Xn+l) 

denote the unique maximal homogeneous ideal in either S or in A. 

Definition 8.15. The socle of A is 

Soc(^) (0 : m) := {/ e A : /m = 0}. 

Definition 8.16. The ring A is Artinian if dime ^ < oo. In that case, we write 

A^C®AiQ)---®Ai, 

with A( 7^ 0. The number £ is the socle degree of A. It is the least number i such 
that m'^+i C /. 

Definition 8.17. The ring A is Gorenstein if it is Artinian and dimcSoc(yl) = 1 
(so Soc(A) = Ai and dime Ag — 1). 

Proposition 8.18. Suppose A is Artinian with socle degree i. Then A is Goren- 
stein if and only if dime Ag — 1 and the pairing given by multiplication 

Ad X Ag^d Ai^C 

is a perfect pairing. 

Proof. See the proof of, and remarks following. Proposition 8.6, |15| . □ 

As an easy corollary, we have 

Corollary 8.19. The Hilbert function of an Artinian Gorenstein ring A is sym- 
metric. That is, if the socle degree of A is I, then 

HAid)=HAii-d) 

for all d. 

Now let S" = C[yi, . . . ,?/„+i], and let S act on S' by treating each Xi as the 
differential operator d / dyi . 

Theorem 8.20. (Macaulay, cf. Theorem 8.1 [15]/' ^'^e "ring A = Sjl is Gorenstein 
with socle degree i if and only if there exists a nonzero g G S*^ such that 

I = ann(g) := {/ G 5 : . . . , d/dyn+i)g = 0}. 

Now consider the case where / is the homogeneous toppling ideal for a sandpile 
graph G with vertices {vi, . . . , w„+i}. Let X — Z{I) be the zero set of / as discussed 
in section 6.1.2 Let a g 5 be a linear polynomial that does not vanish at any point 
of X. For instance, a may be any of the indeterminates, Xi. Restricting the exact 
sequence given by multiplication by a, 

0^ A^ A^ A/{a) 0, 

to each degree d, we find that the Hilbert function for A/ (a) is the first differences 
of the Hilbert function for A, i.e., i?^/(a')(c?) — AHA{d). It then follows from (6.3) 
that A /{a) is Artinian. 
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Definition 8.21. Continuing tlie notation from above, the ring A/{a) is called 
an Artinian reduction of A. Let £ be the socle degree of an Artinian reduction 
of A, and let hd '■— AiJyi(d) for d — 0, . . . , £. Then (/iq: ■ • ■ j hg) is the homogeneous 
h-vector of G (or / or X). 

Remark 8.22. The homogeneous /i-vector and the /i-vector appearing in Defini- 
tion 6.11 are identical in the case the A{vn+i) is in the span of {A{vi) : 1 < i < n} 



(see the discusion after Example 6.12). 

Definition 8.23. We say G is a Gorenstein sandpile graph if its homogeneous 
coordinate ring has a Gorenstein Artinian reduction. We also say that / and X are 
(arithmetically) Gorenstein. 

Remark 8.24. 



(1) Using the notation preceding Definition 8.21 it turns out that if A has 
a Gorenstein Artinian reduction, then every Artinian reduction of A is 
Gorenstein. 

(2) The notion of a Gorenstein ideal is much more general, but requires a dis- 
cussion of the Cohen-Macaulay property, which our toppling ideals (defining 
a finite set of projective points) satisfy automatically (cf. [13]). 

It is well-known that complete intersection ideals are Gorenstein (cf. ^21.8(13)). 
In particular, we have the following. 

Theorem 8.25. Let G be a sandpile graph. If G is a complete intersection, then G 
is Gorenstein. 

Theorem 8.26. Let I be the homogeneous toppling ideal of the sandpile graph G 
having n + 1 vertices. The following are equivalent: 

(1) G is Gorenstein; 

(2) if the minimal free resolution for I is 

then Fn ~ S as an S-module; 

(3) the homogeneous h-vector for G is symmetric. 

Proof. The equivalence of items ([I]) and ([2l is a standard result (cf. [^). The 
equivalence of items ([I]) and ([s]) follows by ]l2j since / is a Cayley-Bacharach ideal 
by Proposition |6.18| □ 



Example 8.27. Let G be as in example ] 7. 4 1 We saw that the last nonzero module 
in the free resolution for I{G)^ is 5'(— (1, 0, 2, 2)), which is isomorphic to S as an 
5'-module. Thus, the caption for Figure [8) stating that G is Gorenstein, is justified 
by ([2| above. 

Define a loopy tree to be a (finite) graph that is formed from a weighted, undi- 
rected tree by adding weighted loops at some (maybe none) of the vertices. 

Theorem 8.28. For an undirected sandpile graph G, the following are equivalent: 

(1) G is a loopy tree; 

(2) G is a complete intersection; 

(3) G is Gorenstein. 
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Proof. Let G = {V, E, s) be a undirected sandpile graph. First suppose that G 
is a loopy tree. Removing any outgoing edges from s leaves a completely wired 
graph having the same homogeneous toppling ideal as G. Hence, G is a complete 



intersection by Theorem 8.13 and hence G is Gorenstein by Theorem 8.25 

We now assume that G is not a loopy tree. Since the lattice ideal of G is 
not affected by loops, for ease of exposition we assume that G has no loops. By 



Theorem 7.7 and Theorem 8.26| ([2]), we have that G is Gorenstein if and only if it 



has a unique minimal recurrent configuration. 

To characterize the minimal recurrent configurations, let ^ be a total ordering 
of the vertices such that for all nonsink vertices v, (i) s ^ v, and (ii) there exists 
u ^ V such that {u, v} G E. Define the configuration by 

c^^v ■= deg(t;) — e V : {u,v} £ E and u -< v}. 

We now invoke Dhar's burnin g algorithm. Let b be the mini mal b urning configura- 



tion for G. By Theorem 2.28 it has script 1, and by Theorem 2.27 a configuration c 
is recurrent if and only if each nonsink vertex fires in the stabilization of 6 + c. Note 
that 6 + c is obtained by starting with c and firing the sink vertex. It follows that 
is a minimal recurrent configuration and that all minimal recurrent configurations 
arise as for some ordering -< satisfying (i) and (ii), above. 

Let C be a (undirected) cycle in G. Choose a path P in G starting at s and going 
to a vertex of G, then traveling around G. To be precise, let ui, . . . ,Ui be distinct 
vertices forming a path in G (so {uf,Ui-^-l} G V for all £) with ui = s and Ui a 
vertex in G. Assume that Ui is the first vertex in the path to be in G. (If s is in G, 
then i — 1.) Next, let Ui, . . . , Uij^j be the vertices in the cycle G, in order. Then P 
is the path ui, . . . , Ui+j. Let be any total ordering satisfying (i) and (ii), above, 
with 

Ul ^1 • • • ^1 U^+j, 

and such that Uk -<i v for all Uk and all vertices v not in P. Let -<2 be any total 
order satisfying (i) and (ii) with 

and such that Uk -<2 v for all and all vertices v not in P. It follows that c^^ 
and distinct minimal recurrent configurations on G. Hence, G is not Goren- 

stein. □ 

By Theorem |5.13[ an Artinian reduction of A for a sandpile graph with absolute 
sink has the set 

{x'^ : c is a superstable configuration of G} 

as a normal basis. It follows that the socle degree ^ of ^ is the maximum of 
the degrees of the superstable configurations of G. Hence, by Theorem 8.26| (|3|, 



sandpile graph with absolute sink is Gorenstein if and only if there exists a bijection 
between the superstable configurations of degree k and those of degree £ — k. 

Lemma 8.29. Let Gi be the graph on a single vertex v and let G2 be a Gorenstein 
sandpile graph. Let G be a wiring of Gi into G2 . Then G is Gorenstein. 

Proof. Let A be the set of superstable configurations on G2 and define the integer 
i :— max{deg(a) : a G A}. Let f : A A he a bijection such that deg(/(a)) = 
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i — deg{a) for all a £ A. Let d := outdeg(u). Since there are no edges from vertices 
of G2 to V in G, the set of superstable configurations on G is 

B {kv + a : a e A and < k < d}. 

Let TO :— max{deg(&) : b E B} ~ £ + d — 1. Define g: B ^ B hy g{kv + a) — 
(c? — 1 — k)v + /(a) where a E A. Then g is a bijection and 

deg{g{kv + a)) = deg((d - I - k)v + f{a)) =£ + d-l- deg{kv + a). 

It follows that G is Gorenstein. □ 



Lemma 8.30. Let Gi be a Gorenstein sandpile graph with absolute sink s and 
let G be a wiring of Gi into the graph on a single vertex v with no edges. Then G 
is Gorenstein. 

Proof. Let A = be the Laplacian matrix for G, and let D be the wiring divisor 
of G. If d is the weight of the edge from s to then 

- -d 

Since \D\ 7^ by the definition of a wiring, there exists some effective divisor 
E D. Thus, we can replace the last column of A^ with 

E 
-d 

without changing the column span, and hence without changing the associated 
lattice ideal. Negating this column and swapping rows, the matrix A^ becomes 

d 



A° := 

which is the restricted Laplacian for a wiring of vertex v into Gi. This graph is 



Gorenstein by Lemma 8.29 □ 



Theorem 8.31. Let Gi and G2 be Gorenstein sandpile graphs with absolute sinks. 
If G is a wiring of Gi into G2, then G is Gorenstein. 

Proof. Let D be the wiring divisor of G. Let G" be the wiring of Gi into the 
graph on a single vertex s, disjoint from the vertices of Gi or of G2, with wiring 
divisor D. Let A' be the set of superstable configurations on G' and define the 



integer £' :— max{deg(c) : c G A'}. Since G" is Gorenstein by Lemma 8.30 there 
exists a bijection f':A'^ A' such that deg(/'(c)) — £' — dcg(c). Let A2 be the 
set of superstables on G2, let £2 = max{deg(c) : c G A2}, and let f2' A2 ^ A2 be 
a bijection such that deg(/2(c)) = £2 — deg(c). 

Clearly, if c is superstable on G, then c\y^ G A2, and c\vi G A'. Conversely, 
if c' G A' and C2 G A2, then the configuration c' + C2 is superstable on G. Let 
A = {cJ + C2 : c' G .4', C2 G A2}, so that A is the set of superstable configurations 
on G, and max{deg(c) : c G .4} = £' + £2 —■ £■ Define the function f : A ^ Ahy 
f{c' + C2) = /'(c') + /2(c2), where c' G A' and C2 G A2. Then / is a bijection, and 
deg(/(c' + C2)) = £ — dcg(c' + C2). Hence, G is Gorenstein. □ 
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